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RESUMEN 

En este articulo se estudia el problema de cargas vivas vehiculares en puentes en Mexico. Se 
discuten los modelos de carga de diseno que han sido utilizados en Mexico y los conceptos de 
impacto y flujo vehicular. Para conocer y entender las solicitaciones que soporta la infraestructura 
mexicana de puentes se realiza el analisis estadistico de los pesos totales y distancias 
caracteristicas de los principales vehiculos de autotransporte provenientes del aforo realizado por 
el IMT. Se estudian dos situaciones extremas: vehiculos en movimiento (con efecto de impacto) y 
vehiculos sin movimiento (sin impacto). Adicionalmente, se evalua el efecto de los vehiculos 
cargados con caracteristicas representativas y los vehiculos de diseno. Se concluye que las 
grandes solicitaciones actuales no pueden representarse con los modelos de carga mas usados y 
que esta conclusion puede hacerse extensiva a cualquier pais que use normas de diseno ajenas a 
sus condiciones locales. 



ABSTRACT 

In this paper the problem of vehicular live vehicular loads on bridges in Mexico is studied. 
Design load models which have been used in Mexico are discussed as are also impact and 
vehicular flow concepts. To find out and understand the demand to which Mexican bridges are 
subjected, a statistical analysis of the total weight and the characteristic distances of heavy 
vehicles is carried out using data provided by the IMT. Two extreme situations are studied: 
vehicles in motion (with impact effect) and vehicles at rest (without impact). Additionally, 
vehicles with typical characteristics and design vehicles are evaluated. It is concluded that the 
current large demands cannot be represented by load models commonly used. This conclusion 
may be applicable to any country that do not take into consideration local conditions. 
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1 INTRODUCTION 



Los vehiculos que conforman el transporte carretero han evolucionado signiflcativamente a lo 
largo de este siglo, presentando caracteri'sticas particulares en cada pais, resultado su 
reglamentacion de los pesos y dimensiones de los vehiculos de carga, la efectividad de su sistema 
de control en el cumplimiento de las normas y sus necesidades para el transporte local. 
Desafortunadamente, no todos los pai'ses poseen informacion actual y confiable sobre las 
caracteristicas del transito y de los vehiculos que circulan por sus redes de caminos. Esta 
informacion es indispensable para la formulacion de un modelo de cargas para puentes que refleje 
las solicitaciones reales a las que estara sometido este tipo de estructuras durante su vida util. Por 
lo que, aun cuando se pueda importar la mayor parte de un codigo de diseno estructural de un pais 
a otro, el modelo de carga de diseno para puentes debe ser definido de manera particular en cada 
pais. 

En Mexico, al igual que en otros paises en desarrollo, no se tiene un codigo propio para el 
diseno estructural de puentes y se han utilizado, para este proposito, las normas estadounidenses 
AASHTO (American Association of State Highway and Transportation Officials). En 
consecuencia, los puentes mexicanos se han disenado con los criterios de analisis y diseno, asi 
como con los modelos de carga viva propuestos por estas normas. El modelo de carga de la 
AASHTO aparecio en la edicion de 1944 y continua sin modificaciones hasta su ultima edicion 
(1996). Desde entonces, el tamano y la capacidad de carga de los vehiculos han evolucionado, 
existiendo incrementos significativos en sus pesos totales que, en general, han sido diferentes 
entre Mexico y Estados Unidos. 

Como resultado de la evolucion de los vehiculos de carga en Mexico, la Secretaria de 
Comunicaciones y Transportes (SCT) publico un reglamento de pesos y dimensiones, en el cual 
se establecieron nuevos pesos y dimensiones permisibles con los que los vehiculos de 
autotransporte pueden transitar en los caminos y puentes de jurisdiccion federal (SCT, 1980). Con 
base en estas normas se han publicado hasta la fecha varias modificaciones a las mismas, entre las 
que resaltan el Proyecto de Norma Mexicana NOM-12-SCT-2-1995 (SCT, 1995) y sus 
modificaciones de 1997 (SCT, 1997). En la norma de 1980 se definieron los vehiculos T3-S3 y 
T3-S2-R4 (apendice A) como los mas pesados; estos vehiculos han sido adoptados, sin ningun 
estudio base, como cargas de diseno por la SCT (Galindo, 1998), sin embargo, no existe ningun 
reglamento donde se definan dichas cargas como cargas vivas de diseno de puentes. 

A partir de 1991, el Institute Mexicano del Transporte (IMT) comenzo a realizar aforos de 
pesos y dimensiones de los vehiculos de carga que circulan por la red mexicana de carreteras. En 
estos estudios se observo que un porcentaje significativo de los vehiculos de autotransporte 
exceden el peso reglamentario. La falta de un reglamento de diseno de puentes en Mexico y la 
disponibilidad de los dates anteriores motivaron la realizacion de la presente investigacion, cuyo 
objetivo principal es el analisis y evaluacion de las cargas vehiculares y sus efectos en los puentes 
mexicanos. 

La falta de control en el cumplimiento de los pesos reglamentarios de los vehiculos de 
autotransporte da lugar a solicitaciones no contempladas por las cargas vivas de diseno actuates. 
Esta falta de control es evidente en el numero significativo de vehiculos sobrecargados que 
circulan en el sistema carretero mexicano. Opiniones dentro del sector transporte en Mexico 
establecen que, con el Tratado de Libre Comercio de Norteamerica (IMT), el peso de los 
camiones mexicanos bajara debido a que en los Estados Unidos los pesos permisibles son 
menores y que los camiones mexicanos que deseen cruzar la frontera norte deberan cumplir las 
normas estadounidenses. Visto objetivamente, lo anterior no necesariamente se cumplira en la 
totalidad de los vehiculos de nuestro pais, ya que solo los transportes transnacionales lo 
cumpliran, mientras que el transporte intemo seguira, sin control, manejando las mismas cargas. 
Esta situacion esta llevando a la infraestructura carretera mexicana a situaciones de deterioro no 
contempladas en el diseno de puentes, al haberse usado cargas no representativas de la realidad 
de este pais. 
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En este articulo se discute el problema de cargas vivas vehiculares en puentes y los modelos 
de diseno que han sido utilizados en Mexico. Se discuten los conceptos de impacto y flujo 
vehicular, los cuales son fundamentales para el analisis de este problema. Se realiza el analisis 
estadistico de los pesos totales y distancias caracterfsticas de los principales vehiculos de 
autotransporte, a partir del aforo realizado por el EMT. Se estudia el problema de cargas vivas 
vehiculares en puentes considerando dos situaciones extremas: vehiculos en movimiento 
(incluyendo el efecto de impacto) y vehiculos sin movimiento (sin impacto). Adicionalmente, se 
evaluan los vehiculos cargados con caracterfsticas representativas y los vehiculos de diseno. Los 
puentes son modelados como vigas simplemente apoyadas y se utilizan los momentos y cortantes 
maximos como parametros de comparacion en el analisis. 



2 MODELOS DE CARGA VIVA EN REGLAMENTOS 

Para poder realizar un analisis critico del problema de cargas vivas vehiculares se requiere, 
primeramente, deftnir dos tipos de cargas vehiculares: 

• Cargas permisibles. Son aquellas que la institucion reguladora del transporte carretero 
establece como los maximos pesos que puede tener cierto vehiculo al circular por los 
caminos de su responsabilidad. 

• Cargas de diseno. Son aquellas que la institucion reguladora de los puentes considera que 
refleja las solicitaciones a las que va estar sometido esta clase de estructuras durante su 
vida util. 

Estos dos tipos de cargas buscan dar solucion a dos problemas completamente distintos, 
aunque relacionados entre si. Una diferencia elemental es que las primeras corresponden a 
vehiculos reales con caracterfsticas admisibles, mientras que las segundas pueden o no ser 
vehiculos reales, por ejemplo los modelos de carga de la AASHTO (1996) no corresponden a 
vehiculos reales sino a modelos representativos de las solicitaciones esperadas. Se mencionan 
estas diferencias ya que la SCT usa las cargas permisibles T3-S3 y T3-S2-R4 como cargas de 
diseno, lo cual, por lo anteriormente expuesto, no es la opcion mas adecuada desde el punto de 
vista conceptual. 

Como se menciono anteriormente, actualmente en Mexico se utilizan las normas AASHTO 
como especificaciones de diseno de puentes con las cargas vehiculares que recomienda la SCT. 
Sin embargo, hay que diferenciar las normas AASHTO 1996, que han sido utilizadas en Mexico 
regularmente, de las normas AASHTO LRFD. Las primeras manejan dos metodos de diseno. 
esfuerzos permisibles y factores de carga. Las segundas estan basadas en el uso de factores de 
carga y resistencia (LRED, por sus siglas en ingles). Estas especificaciones presentan varios 
cambios en relacion con las normas tradicionales, resaltando la utilizacion de factores de carga y 
de resistencia definidos a partir de resultados de estudios estadisticos de la variabilidad de las 
cargas y las resistencias y de la asignacion de niveles de seguridad uniformes (Morcos y Grubb, 
1995). 

A continuacion se describen los modelos de carga viva de las normas para puentes AASEITO 
1996 y AASHTO LRFD 1998, las cargas recomendadas por el Comite de Cargas y Fuerzas en 
Puentes del Comite de Puentes de la Division Estructural de la Sociedad Americana de Ingenieros 
Civiles (ASCE, 1981) y el Reglamento de Pesos y Dimensiones de la SCT. La informacion 
anterior es necesaria para evaluar el problema de cargas vivas vehiculares en Mexico. 

2.1 Normas AASHTO 1996 

La primera carga de diseno empleada en Mexico fue la H-15. Posteriormente, en 1956, se 
utilizaron las cargas HS-15, y a partir de 1972 se emplearon las cargas HS-20 (Fossas, 1991); 
estos modelos de carga estan especificadas en la parte 3.7, “Cargas del Camino de la Seccion 3 
“Cargas”, de las Normas para Puentes de la AASHTO (1996) y fueron por primera vez 
presentadas en su edicion de 1944. 
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Las Normas AASHTO (1996) establecen que las cargas vivas que se deben aplicar sobre la 
calzada de un puente para fines de diseno consisten en las de un camion hipotetico o una carga de 
carril (equivalente a un tren de camiones), definidas ambas para producir acciones maximas sobre 
un claro de un puente. Estas cargas se muestran en las figuras 1 y 2. En la figura 1, la carga H 
corresponde a un camion hipotetico de dos ejes y la carga HS a un tractor de dos ejes con un 
semiremolque de uno. El numero indica el peso bruto de un camion estandar en el caso de la 
carga H, y en el otro, representa el peso bruto del tractor, ambos expresados en toneladas. 
Utilizando la nomenclatura mexicana, definida en el apendice A, el camion H equivale a un C2, 
mientras que el camion HS equivale a un T2-S1. Adicionalmente, las normas AASHTO proponen 
una carga militar compuesta por dos ejes, separados a cuatro pies, con 24,000 libras cada uno. 




1 4.27 m. 

H20-44 3.63 ton 

H15-44 2.72 ton 



14.51 ton 
10.89 ton 



i 



4.27 m 



4.27 m a 9.14 m 



HS20 44 3.63 ton 14.S1 ton 14.51 ton 

HS15-44 2.72 ton 10.89 ton 10.89 ton 



Figura 1. Camiones de diseno, AASHTO (1996) 

8 16 ton para el Momento 
111 79 ton para el Cortante 

♦ 0.95 ton/m 

H20-44 HS20-44 

i <5. 12 ton para el Momento 
8.85 ton para el Cortante 

0.71 ton/m 



H15-44 HS 15-44 



Figura 2. Carga de carril, AASHTO (1996) 



2.2 Normas AASHTO LRFD 1998 

En los Estados Unidos se ha observado que las cargas de diseno HS-20, la carga de lrnea y el 
camion militar, producen momentos y cortantes inferiores a los obtenidos por el transito actual. 
Se ha observado que al pasar los camiones mas pesados sobre sus puentes los momentos y 
cortantes que se obtienen fluctuan entre 1.5 y 1.8 veces los correspondientes a la carga de diseno 
(utilizando la carga HS-20) (Nowak, 1995). Investigaciones recientes (Nowak, 1995) muestran 
que estas diferencias desaparecen si a las cargas del camion de diseno HS-20 o al eje tandem de 
diseno (en el caso de claros cortos) se superpone una carga uniformemente repartida de 9.3 kN/m. 
Con base en estos resultados se llego a la carga de diseno propuesta en las Normas AASHTO 
(1998) (figura 3). 
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14.8 ton 14.8 ton 



3.6 ton 



I i 



.95 ton/m 



1 1 .2 ton 1 1 .2 ton 




.95 tordm 



4.3m 4. 3-9 ,0m 12m 

Figura 3. Camion y carga uniforme, Eje tandem y carga uniforme, AASHTO (1998) 



2.3 Cargas recomendadas por el Comite de Cargas y Fuerzas en Puentes del ASCE (1981) 

Estas recomendaciones fueron resultado de un conjunto de investigaciones encausadas hacia 
el estudio de los modelos de carga viva para puentes de grandes claros. Lo anterior se fundamenta 
en el hecho de que las solicitaciones maximas en un puente de claro corto son funcion del paso de 
un solo vehiculo; sin embargo, para puentes de grandes claros no se cumple esta condition. 
Algunas de las conclusiones de estos estudios fueron (Xanthakos, 1994): 

♦ La carga maxima ocurre cuando los vehiculos no tienen movimiento. 

♦ Esta carga puede ser representada por una carga uniformemente repartida mas una carga 
concentrada. 

♦ La carga promedio disminuye con el incremento de la longitud cargada. 

El Comite propuso que la carga basica de diseno para puentes de grandes claros este definida 
por una carga uniformemente distribuida, U, y una carga concentrada, P. Los valores de U y P 
son funcion de la longitud cargada y del porcentaje de vehiculos pesados. Los valores de U y P 
estan defmidos en ASCE (1981). Dado que estas cargas son las producidas por vehiculos sin 
movimiento, los resultados no se deben afectar por el factor de impacto. 

2.4 Pesos y dimensiones permisibles en Mexico. 

En Mexico se han identificado los vehiculos: C2, C3, T3-S2, T3-S3 y T3-S2-R4 (apendice A), 
como los vehiculos de carga mas comunes que transitan por la red mexicana de carreteras (IMT, 
1995). Entre estos cinco vehiculos, el T3-S3 y T3-S2-R4 (figuras 4 y 5) son los mas pesados, y por 
ello han sido comunmente usados como camiones de diseno para puentes en Mexico. 




+ Dimensiones maximas - Dimensiones minimas 

Se autoriza el transito de estos vehiculos solo por caminos tipo “A ’. 



Figura 4. Vehiculo de carga T3-S3 (SCT, 1980) 
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i- 1.80 m +5.50 m + 1.50 m +5.30 m 

- 3.50 m - 1 .20 m - 4.25 m 

55 ton 9.0 ton 9.0 ten 



+ 1 .50 m + 1 .50 m 

- 1.20 m - 3 .20 m - 1.20 m 

9 JO ton 9J0 ton 9J0ton9J0ton 

PBV775 ton 



+ 5.375 m 
- 4.25 m 



+ 1.50 m+2.80 m 
- 1 .20 m 

9.0 ton 9 JO ton 



+ Dimensiones maximas - Dimensiones mfnimas 

Se autoriza el transito de estos vehi'culos solo por caminos tipo “A”. 



Figura 5. Vehiculo de carga T3-S2-R4 (SCT, 1980) 



La SCT publico en las ediciones del 3 de octubre y 28 de noviembre de 1980 del Diario 
Oficial de la Federacion, el Capftulo XI del Reglamento del Capftulo de Explotacion de Caminos 
de la Ley de Vlas Generales de Comunicacion, que define los valores permisibles de los pesos 
totales y otras caracteristicas de los vehfculos de carga de circulation en Mexico. Esta 
publicacion es la actualization del mismo documento publicado el 24 de diciembre de 1960, 
dando respuesta a la necesidad de poder transportar mayor peso (SCT, 1980). Los valores 
permisibles de los pesos y dimensiones, para distintos tipos de ejes y vehfculos, estan defmidos 
en funcion del tipo de camino (SCT, 1980; SCT, 1995). 

Desde la edition de este reglamento hasta la fecha se han publicado multiples comentarios y 
modificaciones a dichas normas, entre las que se puede resaltar: 

♦ La publicacion, en 1995, del Proyecto de Norma Oficial Mexicana NOM-012-SCT-2-1995 
(SCT, 1995). 

♦ La publicacion del 7 de enero de 1997 presenta algunas modificaciones a la Norma Oficial 
Mexicana NOM-012-SCT-2-1995, (SCT, 1997). 

Los pesos permisibles correspondientes a las ediciones reglamentarias de 1980 y 1995 se 
resumen en la tabla 4. 



3 CONDICIONES DE CARGA EN EL PUENTE INDUCIDAS POR EL FLUJO 
VEHICULAR 

Para el analisis estructural de un puente se debe tomar en cuenta que la superestructura va a 
estar sometida a varios tipos y combinaciones de carga. En la definition de las cargas se debe 
considerar la posibilidad de que los carriles (dos o mas) esten ocupados simultaneamente, y que 
ademas existan posibles variaciones en la velocidad de los vehfculos que transitan sobre el 
puente. 

Este ultimo punto resulta de interes analizarlo, para ello se puede mencionar que existen dos 
situaciones extremas: una es cuando el vehiculo esta transitando a la velocidad de proyecto o de 
uso, y otra cuando la velocidad es nula, por ejemplo cuando el transito en el puente esta 
congestionado. Mientras mayor sea la velocidad, la distancia entre vehfculos sera mayor, 
consecuentemente la densidad vehicular (numero de vehfculos en una longitud determinada) sera 
menor (Xanthakos, 1994). Se sabe tambien que al aumentar la velocidad de los vehfculos se 
incrementa el efecto del impacto (efecto dinamico de la carga vehicular), por lo que un vehiculo 
con velocidad nula o casi nula no produce impacto. 
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A partir de lo anterior se puede definir un criterio de analisis para el estudio de las cargas 
vivas en un carril de un puente. Este criterio consiste en encontrar la situacion mas critica al 
considerar las siguientes dos situaciones extremas: 

♦ La ocurrencia de uno o mas vehiculos recorriendo el puente con una velocidad modal 
(representativa), incluyendo el efecto de impacto. 

♦ La ocurrencia de dos o mas vehiculos con velocidad nula o casi nula sobre el puente, sin tomar 
en cuenta el efecto de impacto y teniendo una distancia entre vehiculos menor a la considerada 
en el primer caso. 



3.1 Impacto 

El factor de impacto se define como la relacion de la respuesta dinamica maxima entre la 
respuesta estatica maxima menos uno. La modelacion del impacto es un problema complejo, ya 
que participan multiples variables entre las que se pueden citar (figura 6): 

♦ Propiedades dinamicas del vehiculo, las cuales dependen de su masa, amortiguamiento, 
dimensiones, etc. 

♦ Propiedades dinamicas del puente, que dependen de su masa, rigidez a flexion, claro, etc. 

♦ Velocidad de recorrido del vehiculo. 

♦ Rugosidad del pavimento, la cual esta determinada por las condiciones de la superficie de 
rodamiento, tanto en el puente como en sus accesos. 

El Codigo de Ontario en 1983 (Bakht y Jaeger, 1987) especifica la carga dinamica como una 
funcion de la frecuencia natural de vibracion del puente. Esta especificacion, al incluir las 
propiedades dinamicas del puente, se toma conceptualmente mas completa, aunque mas 
compleja, en comparacion con las correspondientes a otros codigos, como es el de las Normas 
AASHTO (1996), en donde se considera que el impacto es funcion solamente del claro del 
elemento cargado. Estas especificaciones definen, en el apartado 3.8.2, la formula de impacto 
como: 



Z, + 50 



en donde: 

/= es la funcion del impacto (maximo 30%) 

L= es la longitud de la parte del claro cargado en pies 




Propiedades dinamicas del puente: 
Rigidez, masa, amortiguamiento. 



Figura 6. Variables que intervienen en la determination del impacto 



La modelacion del impacto es dificil de lograr por medio de expresiones sencillas que reflejen 
resultados acordes con la realidad. Una tendencia reciente ha sido el uso de mediciones de campo 
para establecer las solicitaciones reales causadas por dicho efecto. Este proceder se sugiere en la 
primera edition de las Normas LRFD de la AASHTO (1998), en donde se establecen valores de 
impacto independientemente de las caracterfsticas del puente: 
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Tabla 1. Valores de “Impacto”, Normas LRFD, AASHTO, 1998 



Componente 


Impacto 


Juntas de la cubierta. Todos los estados llmite 


75% 


Otros componentes: 

* Estados llmite de fatiga y fractura 

* Todos los demas estados llmite 


15% 

33% 



Desde el punto de vista de los autores, se debe tomar con cuidado el resultado de las 
investigaciones realizadas sobre este tema en otros pafses, ya que cada pals posee condiciones 
particulares; por ejemplo en Mexico: la masa de los vehlculos es grande (sobrepeso), la 
suspension de los vehlculos de carga presentan gran variedad en la rigidez y el amortiguamiento 
(desgaste por el sobrepeso, menor mantenimiento) y las distintas condiciones de la superficie de 
rodamiento (pavimentos y terraplenes en malas condiciones, asentamientos diferenciales a lo 
largo del puente, etc.). 

3.2 Flujo Vehicular 

La ingenierla de transito ha desarrollado diversos modelos, macro y microscopicos, para el 
analisis del flujo vehicular. Estos modelos requieren de la definicion de tres variables principales 
(Papacostes y Prevedouros, 1993; Cal y Mayor y Cardenas, 1994): 

1) Intensidad (q): es la frecuencia con la cual pasan los vehlculos por un punto o seccion 
transversal de un carril o calzada. En otras palabras, la intensidad es el numero de vehlculos N 
que pasa en un tiempo T, (q=N/T). La intensidad es un parametro que varfa a lo largo del dla. 

2) Densidad vehicular (k): es el numero de vehlculos N que ocupan una longitud especlfica d 
(k=N/d). 

3) Velocidad (v): es la distancia que recorre el vehlculo en un tiempo determinado (v=d/T). 

Las tres variables anteriores estan relacionadas entre si mediante la ecuacion fundamental del 
flujo vehicular, cuya forma general es: 

q = vk 



Uno de los modelos macroscopicos del flujo vehicular es el lineal. Este modelo es sencillo y 
de facil interpretacion, ya que establece la existencia de una relacion lineal entre la densidad 
vehicular y la velocidad. La densidad disminuye a medida que aumenta la velocidad y viceversa 
ya que son inversamente proporcionales. 

La Ingenierla de Transito ha recurrido, frecuentemente, a la utilizacion del proceso de Poisson 
para modelar el flujo vehicular (libre) en un punto determinado (Benjamin y Cornell. 1970). La 
distribucion de probabilidades de una variable aleatoria x asociada a un proceso de Poisson que 
representa el numero de eventos que ocurre en un intervalo de tiempo t, esta dada por: 

, ^ e~\Aty 

p(x,At) 

x! 

en donde At es el promedio de eventos en un intervalo dado. 

Para el analisis del flujo vehicular se evaliia la probabilidad de que x vehlculos sean 
observados en un intervalo de tiempo t, dada una intensidad de transito At (vehlculos/tiempo) 
(Benjamin y Cornell, 1970). 

La distribucion de Poisson, que ha sido utilizada para el analisis de problemas de transito. no 
describe adecuadamente la distancia entre vehlculos. Se han realizado varias propuestas entre las 
que destaca la de Andre Schuhl (Kell, 1962), quien supone que el flujo vehicular esta divido en 
dos grupos: 

♦ Vehlculos que fluyen libremente. 

♦ Vehlculos restringidos (vehlculos influenciados por el vehlculo que circula adelante). 
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Cada uno de estos grupos tiene distintas medias y obedecen a alguna distribucion de Poisson. 
La distribucion total esta compuesta por la suma de las dos subdistribuciones y se establece por 
medio de la ecuacion: 

t-A t-t 

p(h>t) = (\-a)e T '~ A + (a)e Tl ~ r 

en donde: 

p(h>t) es la probabilidad de que el tiempo entre vehiculos (h) sea mayor o igual que el 
tiempo (/) 

a es proporcion de vehiculos de circulacion restringida en el flujo vehicular 
(1 - a) es la proporcion de vehiculos de circulacion libre en el flujo vehicular 
j es el tiempo promedio entre los vehiculos de circulacion libre 
j es el tiempo promedio entre los vehiculos de circulacion restringida 
A es el tiempo mi'nimo promedio entre los vehiculos de circulacion libre 
r es el tiempo mi'nimo promedio entre los vehiculos de circulacion restringida 
e es la base de los logaritmos naturales 

En la revision microscopica del flujo vehicular se analiza el comportamiento de dos vehiculos, 
uno tras otro, que se encuentran en un flujo sin interrupciones extemas (por ejemplo, 
entronques). Como una regia general, debe existir una separacion entre vehiculos tal que si existe 
una desaceleracion repentina por parte del vehiculo delantero, el que le sigue debe tener una 
distancia suficientemente amplia para detenerse sin chocar con dicho vehiculo, figura 7 
(Papacostas y Preveduoros, 1993). 
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V : velocidad irucial de los dos vehiculos 
di : desaceleracion del vehiculo 1 
d 2 : desaceleracion del vehiculo 2 
6 : tiempo de percepcion-reaccion del 

conductor del vehiculo 2 
X 0 : distancia de segundad cuando los 
vehiculos estan detemdos 
L : longitud del vehiculo 



Figura 7. Frenado de dos vehiculos consecutivos 



En esta figura la distancia de frenado del vehiculo 2 esta conformada por la distancia de 
percepcion-reaccion y la de frenado: 

V 2 



X, = VS + - 



2 d. 



La distancia de percepcion-reaccion esta defmida por la velocidad del vehiculo y el tiempo 
percepcion-reaccion. El tiempo percepcion-reaccion esta formado por cuatro acciones: 
percepcion, identification, decision y reaccion (McShane y Roess, 1990). Del estudio dirigido 
por Johannseon y Rumor, para evaluar el tiempo de respuesta ante frenado sorpresivo se obtuvo 
que el tiempo de reaccion fluctua entre 0.3 y 2.0 segundos, con un valor medio de 0.66 segundos 
(McShane y Roess, 1990; Papacostas y Preveduoros, 1993). Tradicionalmente, en el diseno 
geometrico de caminos se maneja un tiempo de reaccion de 2.5 segundos. 

Considerando que los dos vehiculos tienen la misma capacidad de frenado (frenan en la misma 
distancia), la distancia entre vehiculos es: 



S - L = VS + X„ 
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4 PROCEDIMIENTO DE ANALISIS 



Los puentes pertenecen a la clase de estructuras que estan sujetas a cargas dinamicas, es decir, 
la posicion de las cargas varfa con el tiempo. En este tipo de estructuras las solicitaciones 
intemas varian a causa de la magnitud y posicion de las cargas extemas; por lo tanto, para disenar 
los elementos estructurales se debe determinar la posicion de las cargas que produce las maximas 
solicitaciones 

En este trabajo se evaluan dos tipos de elementos mecanicos maximos: (1) momentos 
flexionantes y (2) fuerzas cortantes. Estos dos parametros son los puntos de comparacion en los 
analisis realizados. Para obtener los valores maximos de estos elementos mecanicos producidos 
en una viga simplemente apoyada, modelo aplicable a los claros de la mayorfa de los puentes, se 
elaboro un programa de computadora cuyas hipotesis y diagrama de flujo se presentan en 
Fernandez (1997). 

Para el analisis de los dos escenarios de carga, vehfculos en y sin movimiento, se requiere 
determinar la distancia entre vehfculos. Para el caso en el que la velocidad es aproximadamente 
nula se propone, con base en observaciones de campo, una separacion de 4m. entre el eje trasero 
del primer vehfculo y el eje delantero del segundo vehfculo. Para el otro escenario se recurre a la 
distribucion propuesta por Schuhl (Kell, 1962), para lo cual se selecciona los siguientes 
parametros: 

L=600 vehfculos/hora (trafico moderadamente intenso) 

A=0.9 

t= 1.2 

<2=0.3 (vehfculos restringidos) 
l-a= 0.7 (vehfculos no restringidos) 

AT, =6.90 

K 2 = 1.94 

Dada una probabilidad del 95% de que la distancia entre vehfculos sea mayor o igual a un 
tiempo <, por tanteos se obtiene que: 

p(h>t)=95% para <=1.28 seg. (redondeando <=1.3 seg.) 

Considerando una velocidad del transito (moderada) de 80 km/hr (22.2 m/seg.), la distancia 
entre vehfculos es: 

d = vt = (22.2 m/seg.)(1.3 seg.) = 28.89 m. 

Aumentando la distancia debido al vuelo de los vehfculos, la distancia entre ejes resulta: 
d = 3 1 m. 

El tiempo que define la distancia entre vehfculos es 1.30 seg., que corresponde al 52% del 
tiempo de percepcion-reaccion de seguridad de 2.5 seg. Este tiempo es mayor en 97 % a la media 
de 0.66 seg. obtenida del estudio percepcion-reaccion de Johannseon y Rumor (Papacostas y 
Prevedouros, 1993). Lo anterior muestra que el resultado obtenido por la formula de Schuhl 
corresponde a un valor intermedio entre el tiempo medio percepcion-reaccion (distancia minima) 
y el tiempo de seguridad. 

Con el fin de representar un escenario real del transito se usa la base de datos del aforo de 
1993. por considerarse la mas confiable de las existentes, la cual fue revisada y depurada para el 
analisis de los principales vehfculos. Esta base de datos se ordeno por estacion y tiempo de 
arribo, para tener la secuencia real del paso de vehfculos por el camino. 

Entre los camiones de carga circulan, simultaneamente, automoviles y camionetas. De los 
aforos realizados en la red mexicana de carreteras, se observa que alrededor del 30 % de los 
vehfculos son de carga (IMT, 1997). Esta situacion se puede establecer mediante la generacion de 
"umeros aleatorios, de manera que se introduzca en la base de datos la presencia de los 
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automoviles. El paso de vehfculos por un punto se reproduce por medio de la generacion de 
numeros aleatorios utilizando una distribucion de Bernoulli. En el caso de ser 1 el numero 
aleatorio, se incluye un camion procedente de la base de datos. En el caso contrario, es decir 0, se 
incluye un vehiculo pequeno, el cual fue definido con dos ejes de 750 kg y una distancia entre 
ejes de 250 cm, dando lugar a un vehiculo de 1.5 ton de peso total. 

Para representar los dos escenarios se realizaron cuatro analisis por cada estacion, como se 
describe a continuacion: 

♦ Analisis 1. Simulacion de un vehiculo en movimiento. 

♦ Analisis 2. Simulacion de dos vehfculos en movimiento. 

♦ Analisis 3. Simulacion de dos vehfculos con velocidad practicamente nula. 

♦ Analisis 4. Simulacion de tres vehfculos con velocidad practicamente nula. 

Es importante realizar los cuatro analisis debido a que no se conoce, a priori , que condicion 
de carga dominara en claros de mediana longitud. 

5 ANALISIS ESTADISTICO DE LOS DATOS 

En Mexico se han llevado a cabo estudios de aforo vehicular previos a los realizados por el 
IMT; entre los estudios realizados se tiene la investigacion efectuada por el Instituto de Ingenierfa 
de la UNAM (Magallanes y Negroe, 1977). En esta investigacion se establecio una estacion 
muestreadora o pesadora dinamica de transito que, sin interferir con la operacion normal de la 
carretera, permitio clasificar los vehfculos segiin su tipo, velocidad, distancia entre ejes y peso 
por eje. Entre las conclusiones de este trabajo se presentan los siguientes puntos relevantes: 

♦ Aproximadamente 40% de los ejes tandem y de 15 a ”5% de los ejes sencillos traseros, llevaban 
cargas que exceden el lfmite legal. 

♦ Aproximadamente 40% de los vehfculos T3-S2 llevaban cargas superiores a la legal. De los 
demas, menos de 10% iban sobrecargados. Por no existir entonces, en este estudio no se 
analizaron los camiones T3-S3 y los T3-S2-R4. 

5.1 Estudio de pesos y dimensiones del Instituto Mexicano del Transporte 

El Instituto Mexicano del Transporte realiza desde 1991 un proyecto para evaluar los pesos y 
dimensiones de los vehfculos que circulan en la red mexicana de caminos. Para este proposito se 
han instalando estaciones de aforo en las principales carreteras. El procedimiento de aforo 
utilizado consiste en colocar senalamientos a 500m. de la estacion y 100m antes se encausa el 
trafico para encuestarlo. Posteriormente, una brigada de trabajadores mide las dimensiones de los 
vehfculos, toma datos sobre algunas de sus caracterfsticas y encuesta al chofer sobre la carga 
transportada. Por ultimo, 50 m. despues de la brigada se localiza una pesadora dinamica donde se 
registra la fecha, hora, velocidad del vehiculo, peso de los ejes y distancias entre ellos (IMT, 
1992). 

A continuacion se presentan tablas con algunos resultados estadfsticos interesantes, producto 
de los aforos realizados en 1991, 1992 y 1993. 

La tabla 2 presenta la composicion del transito de vehfculos de carga que circulan en la red 
mexicana. Resalta el hecho de que esta composicion esta comprendida en su mayorfa (98%) por 
vehfculos: C2, C3, T3-S2, T3-S3 YT3-S2-R4. 

En la tabla 3 se presenta el peso bruto vehicular promedio de los vehfculos vacfos. Este peso es 
tambien conocido como “p eso de la tara”. La tabla solo presenta los pesos medios correspondientes 
a los vehfculos mas representativos. 

Es importante aclarar que estas estadfsticas son resultado de aforos realizados en diferentes 
lugares; por lo tanto, la composicion vehicular y el volumen de transito medidos son diferentes. 
De igual modo, la informacion presentada en estos informes no permite conocer la evolucion de 
las caracterfsticas de los vehfculos de carga a traves del tiempo. Finalmente, el numero de 
estaciones aforadas fue distinto en cada afio, siendo 10 en 1991, 3 en 1992 y 15 en 1993. 
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Tabla 2. Porcentaje por tipo de vehlculo de carga (IMT, 1992; IMT, 1995) 



TIPO 


1991 


1992 


1993 


PROMEDIO 


C2 


35 


38.8 


49.4 


41.8 


C3 


22 


20.0 


18.1 


20.0 


C4 


1 


0.4 


0.2 




T2-S1 


0 


0.1 


0.1 




T2-S2 


1 


0.2 


0.2 




T2-S3 


1 


0.2 


0 




T2-S1-R2 


0 


0 


0 




T3-S2 


24 


16.3 


15.3 


18.5 


T3-S3 


15 


21.1 


14.3 


16.8 


T3-S1-R2 


0 


0.1 


0 




T3-S2-R2 


0 


0.1 


0.1 




T3-S2-R3 


0 


0.1 


0.1 




T3-S2-R4 


2 


1.7 


1.5 


1.7 


C3-R2 


0 


0 


0 




C3-R3 


0 


0.1 


0 




OTROS 


0 


0.8 


0.5 




TOTAL 


101* 


100 


100 





* El redondeo de la informacion hizo que el total no sume cero. 



Tabla 3. Peso bruto medio de los vehlculos vaclos 



TIPO 


PESO 


C2 


4 ton. 


C3 


9 ton. 


T3-S2 


17 ton. 


T3-S3 


19 ton. 


T3-S2-R4 


27 ton. 



5.2 Analisis de la informacion registrada en 1993 

La base de datos se obtuvo en formato de hoja de calculo (Excel). Los datos se organizaron en 
funcion del tipo de camion, utilizando el numero de ejes y la distancia entre ellos para agrupar 
correctamente estos datos. Se eligio trabajar con los camiones C2, C3, T3-S2, T3-S3 y T3-S2-R4 
por ser los vehlculos que tienen los mayores porcentajes en la composicion del transito de carga 
en Mexico. 

5.2.1 Peso Total. 

En la figura 8 se presentan los histogramas de peso total de los cinco principales vehlculos de 
autotransporte, y en la tabla 4 se muestra un resumen del analisis estadlstico de estos 
histogramas. 
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Figura 8. Histogramas de pesos totales 
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Tabla 4. Resumen del analisis estadfstico de los pesos totales. 



VEHICULO 


PESO 


PERMISIBL 

E 


MEDIA 


DESVIACIO 

N 


NUMERO 

DE 


MODA 




1980 


1995 




ESTANDAR 


DATOS 




C2 


15.5 


17.5 


7.866 


4.786 


72,987 


3.5, 13.5 


C3 


23.5 


26.0 


17.216 


8.663 


22,273 


9.5, 16.5, 26.5 


T3-S2 


41.5 


44.0 


30.244 


13.627 


15,941 


16.5, 33.5, 

40.5 


T3-S3 


46.0 


48.5 


42.666 


21.684 


15,391 


18.5, 38.5, 

59.5 


T3-S3-R4 


77.5 


66.5(72.5*) 


58.683 


28.183 


1,530 


26.5,78.5 



*Peso permisible cumpliendose ciertas condiciones (revisar referencias SCT, 1995 y 1997) 



5.2.2 Distancia Caracteristica 

En esta seccion se analiza la distancia caracteristica, que es definida como la distancia entre 
ejes que mas varia en un vehiculo. Por ejemplo, en el caso del vehiculo C3 existen dos distancias 
entre ejes; la primera (distancia 1-2) es variable, mientras que la segunda corresponde a la 
distancia entre dos llantas del mismo eje (distancia 2-3), la cual es menos variable que la anterior, 
por lo que la distancia 1-2 es la distancia caracteristica. En la figura 9 se presentan las 
definiciones de las distancias caracteristicas para los vehiculos analizados. 
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Figura 9. Distancias caracteristicas 
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En la figura 10 se presentan los histogramas correspondientes a las distancias caracteristicas 
de los cinco vehiculos estudiados. En la tabla 5 se presenta un resumen los datos estadlsticos 
correspondientes. 



C2 



C3 




DISTANCIA CARACTERIST1CA (CM) 
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T3-S2 



T3-S3 




T3-S2-R4 



T3-S2-R4 



n n 
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Figura 10. Histogramas de las distancias caracteristicas 
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Tabla 5. Resumen del analisis estadfstico de las distancias caracterlsticas. 



Vehfculo 


Permisible 


Media 


Desviacion Estandar 




Moda 






1980 






1 


2 


3 


C2 


310-785 


452.5 


124.4 


337.5 


487 5 


487.5 


C3 


425-785 


512.7 


79.2 


487.5 


612 5 




T3-S2 


425-930 


819.7 


144.2 


412.5 


8125 




T3-S3 


425-805 


653.5 


106.9 


387.5 


687.5 




T3-S2-R4 


425-530 


648.7 


144.8 


487.5 


737.5 






425-537.5 


571.1 


121.2 


462.5 


637.5 





5.3 Analisis de los resultados obtenidos 

Los histogramas de pesos totales muestran como todos los vehfculos de carga analizados 
presentan un gran numero de unidades excedidas en peso, salvo el C2 que presenta solo un 
pequeno porcentaje de vehfculos con carga mayor a la permisible. En contraste, en los 
histogramas de distancia caracterfstica se observa que unicamente los vehfculos T3-S2 y T3-S2- 
R4 presentan distancias caracterfsticas superiores a las permisibles. 

Hwang y Nowak (1991) mencionan que la distribucion de los pesos en los trailers de cinco 
ejes es bimodal. Una distribucion es bimodal cuando se presentan dos intervalos modales (de 
mayor frecuencia), es decir, dos jorobas en el histograma. La primera joroba corresponde a los 
vehfculos vacfos y la segunda se refiere a los vehfculos cargados. Por otro parte, Imbsen (1992) 
senala que el histograma se puede subdividir en cuatro (con sus correspondientes frecuencias 
modales): vehfculos vacfos, parcialmente cargados, completamente cargados y sobrecargados. Es 
importante resaltar que la realizacion de una subdivision objetiva del diagrama es diffcil, debido a 
que la agrupacion se lleva a cabo por tipo de vehfculo a partir del numero de ejes y la distancia 
entre los mismos, pero cada tipo de vehfculo presenta caracterfsticas muy variables entre si, por 
ejemplo: la potencia del motor, la capacidad del chasis, la capacidad de las llantas, etc. 



6 EVALUACION DE LOS EFECTOS DE LAS CARGAS VEHICULARES EN LOS 
PUENTES 

Antes de realizar la evaluacion de las cargas vehiculares reales, se evaluaron las solicitaciones 
producidas en una viga simplemente apoyada por los vehfculos de diseno que han sido utilizados 
en Mexico. En la figura 1 1 se presentan las graficas de los momentos y cortantes (normalizados 
con respecto a la carga HS-20) producido por las cargas HS-15, HS-20. T3-S3 y T3-S2-R4. Es 
importante mencionar que las cargas HS estan conformadas por un camion (que domina claros 
cortos) y una carga de lfnea (que domina claros largos), por lo que se observa un cambio de signo 
en la pendiente de los vehfculos T3-S3 y T3-S2-R4. Los vehfculos T3-S3 y T3-S2-R4 
corresponden a los vehfculos de carga reglamentarios (1980) con pesos maximos y distancias 
mfnimas permisibles para caminos tipo A (condicion crftica). 




Figura 11. Solicitaciones producidas por los vehfculos de diseno 
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De estas figuras se puede concluir: 

1) En el claro de 45 y 40m se obtienen los maximos momentos y cortantes producidos por 
los camiones T3-S3 y T3-S2-R4, llegando a alcanzar relaciones de 1.32 y 2.00 para el 
momento, asi como 1.31 y 1.94 para el cortante. 

2) Los efectos de la carga de los camiones T3-S3 y T3-S2-R4 no se pueden representar por 
medio de cargas tipo “HS", a menos que se discretice la grafica y se propongan varias 
cargas “ HS ” que varien como funcion del claro. 

3) Los puentes con claros menores a 12m presentan valores razonablemente similares (hasta 
17%) entre los camiones T3-S3, T3-S2-R4 y HS-20. 

Adicionalmente se puede observar que el modelo de carga T3-S2-R4 es poco practico para ser 
empleado como carga de diseno, ya que esta conformada por 9 cargas puntuales (ejes). 

En la figura 12 se evaluan los efectos producidos por los automoviles (con las caracteristicas 
utilizadas en esta investigacion) y los camiones C2, C3, T3-S2, T3-S3 y T3-S2-R4 con 
caracteristicas modales relacionadas con grandes pesos (tabla 4). Ambas graficas estan 
normalizadas respecto a la carga HS-20. Los pesos y distancias entre ejes se resumen en la tabla 
6. Estas graficas son utiles para comprender los elementos mecanicos promedio que genera cada 
tipo de vehiculo cargado. 

Tabla 6. Caracteristicas modales de los vehiculos cargados 



Vehiculo 


Peso por eje (1, N) 
(kg) 


Distancia entre ejes 
(1-2, N-1=>N) 
(cm) 


Peso 

total (kg) 


Automovil 


750, 750 


250 


1,500 


C2 


4,500 - 9,000 


587.5 


13,500 


C3 


5,000 - 10,750 - 10,750 


612.5-125 


26,500 


T3-S2 


5,000 - 8,875 - 8,875 - 8,875 


450- 130-812.5- 120 


40,500 


T3-S3 


5,000 - 12,100 - 12,100 - 
10,100- 10,100- 10,100 


450- 130-687.5 - 120- 120 


59,500 


T3-S2-R4 


5,000 - 9,187 - 9187 - 9,187 - 
9187 -9,187 -9187 -9,187-9187 


450 - 130 - 737.5 - 120 - 270 
- 120-637.5- 120 


78,500 



Nota: La numeracion de los ejes es de la cabina del conductor hacia la parte posterior del 
vehiculo. 




CLARO DEL PUENTE (m) CLARO DEL PUENTE (m) 



Figura 12. Momentos y cortantes generados por vehiculos modales 
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Las graficas de estas figuras permiten concluir: 

1. Se presentan elementos mecanicos superiores a los contemplados en las normas para claros entre 
6 my 12 m generados por el T3-S3, cuyo peso excede al permisible. 

2. En claros menores a 8 m se observan momentos elevados que estan definidos por el peso de un 
eje. 

3. Mientras que entre los vehfculos de diseno domina el T3-S2-R4 para casi todos los claros 
analizados, entre los veln'culos modales domina el T3-S3 para claros menores a 25m y para los 
mayores a 25m domina el T3-S2-R4. 

4. Los claros de 14 y 16m son los que presentan menores momentos relativos en comparacion con 
la carga HS-20, siendo del orden de 1.22. Los cortantes relativos son menores en el claro de 10m. 
siendo del orden de 1. 1 7. 

5. Los claros entre 40 y 45m presentan relaciones de cortantes y momentos mayores en 
comparacion con la carga HS-20, siendo de 1.82 y 1.88 respectivamente. 

A continuacion se presentan los histogramas correspondientes al Analisis 1 para un claro de 
14 metros. A estas graficas se les quitaron los automoviles v se dejaron solamente los cinco 
principales vehiculos de carga. Las li'neas verticales corresponden a: la media, la media mas una 
desviacion estandar (D.S.), la media mas dos D.S. y la media mas tres D.S. Los histogramas 
correspondientes a otros claros son similares en forma (Fernandez, 1997). 




MOMENTO (KG CM) 



[ M| FRECUENCIA -»-FREC. ACUM. (%) j 




Figura 13. Histogramas de momentos y cortantes para 14m. de claro. Analisis 1 



En la figura 14 se grafica la media mas n (n=0, 1, 2 o 3) desviaciones estandar del analisis 1 
para varios claros, desde 6 hasta 70 m. Estos valores son para una viga simplemente apoyada y 
estan normalizados respecto a la carga HS-20. El valor del coeficiente de variacion varfa de 0.78 
para claros de 6m. hasta 0.93 (aproximadamente) para claros de 70m. 




Figura 14. Analisis estadfstico de la Analisis 1 para momentos y cortantes 



FREC. ACUM. (% ) 
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En la presente investigacion se compararon los efectos generados por los dos escenarios 
extremos, utilizando una media mas tres desviaciones estandar de los momentos y cortantes 
ohtenidos: 

♦ Caso 1. Vehiculos en movimiento. Analisis 1 y 2. 

♦ Caso 2. Vehiculos sin movimiento. Analisis 3 y 4. 

Cuando se considero el efecto de impacto de las normas AASHTO 1996, el cual disminuye al 
aumentar el claro, los valores del caso 1 dominaron hasta los 45m. aproximadamente y para 
claros mayores el caso 2 fue ligeramente superior. Cuando se considero el efecto de impacto de 
las normas AASHTO LRFD 1998, i.e. constante para todos los claros, el caso 1 domino todo el 
rango estudiado (6-70m.). 

Consecuentemente, el resultado del analisis de los dos escenarios de carga depende 
principalmente de la seleccion del factor de impacto. El analisis 3 no domino en ninguno de los 
casos. De los resultados aqui presentados se puede concluir que el analisis de un puente, dentro 
del rango de 6 a 50 m., se puede llevar a cabo mediante el paso de un solo vehiculo critico sin 
caer en errores significativos. 

Finalmente, en la figura 15 se presenta la comparacion de los analisis realizados cn este 
trabajo. En estas figuras se graftca la carga HS-20, la envolvente de las cargas de diseno o 
permisibles T3-S3 y T3-S2-R4 (SCT, 1980), la envolvente de los vehiculos con caracteristicas 
modales y la media mas tres desviaciones estandar del caso 1 (Analisis 1 y 2). 




— a — Media ♦ 3 Deav S 
— «k — Vehiculos modales 
— — HS-20 

— *t— Vehiculo perrmibte 




Media + 3 Desv S 
[—■* — Vehiculos modales 
HS-20 

Vehiculo permislble 



CLARO DEL PUENTE (m) 

Figura 15. Resumen de los casos estudiados 



En la figura 15 se observa que la forma de la curva de vehiculos cargados con caracteristicas 
modales (solicitaciones reales) difiere de la curva de vehiculos permisibles (modelo de diseno). 
Es importante mencionar que la forma de la curva de solicitaciones maximas pudiese diferir de la 
curva de vehiculos modales, sin embargo, las solicitaciones maximas presentaran grandes 
diferencias respecto a los modelos actuales de diseno. El desarrollo de modelos a partir de 
resultados estadisticos, unicamente, (i.e. una media mas cierto numero de veces la desviacion 
estandar) es inadecuado. La afirmacion anterior se realiza con base, entre otros motivos, en que 
tanto los coeficientes de variacion de la distribucion de solicitaciones como el tipo o parte de 
camion o camiones que domina un claro especifico, varia con respecto a cada longitud de claro 
de puente, aunque los histogramas presenten globalmente formas similares para distintos claros. 



7 CONCLUSIONES 

En esta investigacion se analiza y evalua el problema de cargas vivas vehiculares y sus efectos 
en puentes en Mexico. Esto se lleva a cabo utilizando la base de datos del aforo vehicular de 
1993 realizado por el IMT. El problema es abordado desde sus fundamentos a traves de dos 
situaciones extremas: vehiculos en y sin movimiento. Lo anterior nos condujo a revisar los 
conceptos de flujo vehicular e impacto. De los resultados de esta investigacion se puede concluir: 
1- Como parte del desarrollo de unas normas para el diseno de puentes en Mexico, se requiere 
establecer una carga vehicular de diseno que cubra las solicitaciones reales dadas por los 
vehiculos que circulan en la red nacional de carreteras. No es aceptable establecer una carga 
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vehicular de diseno basada en cargas vehiculares permisibles, ya que estas presentan variaciones 
o ajustes que responden a causas diferentes a la seguridad estructural de los puentes. 
Adicionalmente, es deseable contar con una carga de diseno sencilla y no una complicada como 
la T3-S2-R4. 

2 - Un gran porcentaje de los vehiculos cargados de autotransporte presenta pesos totales que 
exceden a los permisibles, lo que muestra una falta de control en los pesos vehiculares de nuestra 
red de carreteras, quedando como limite de cargas el criterio del usuario y la capacidad del 
vehfculo. Estas solicitaciones no estan contempladas por las cargas vivas de diseno actuales. 

3. - Para claros simplemente apoyados mayores a 8 m, la carga permisible que domina es la del 

camion T3-S2-R4 (SCT, 1980). La envolvente de las cargas permisibles no puede ser 
representada por medio de una carga tipo “HS”, a menos que se discretice la grafica y se 
propongan varias cargas “HS”, las cuales varfen en funcion del claro. 

4. - Del analisis de los vehiculos cargados con caracterfsticas modales, se obtuvo que para los claros 

menores a 25 m domina el camion T3-S3, y para mayores a 25 m rige el T3-S2-R4, lo cual 
difiere del modelo de diseno (permisible). El T3-S2-R4 permisible genera valores superiores al 
modal debido a que su distancia caracterlstica es menor a la modal. 

5. - Los puentes de claros cortos pueden presentar problemas de resistencia y/o fatiga (puentes de 

acero), ya que todos los modelos de carga hasta ahora utilizados subestiman las solicitaciones. 
En puentes de concreto presforzado se puede presentar agrietamiento debido a solicitaciones 
superiores a las contempladas, que aumenta el fenomeno de corrosion en la estructura. 

6. - Para efecto de estudios futuros es deseable contar con mediciones realizadas en distintos tipos de 

caminos y en los mismos caminos realizar mediciones periodicas para conocer la evolucion de 
los vehiculos de autotransporte y los pesos de los mismos. Adicionalmente, es necesario realizar 
mediciones sin interrumpir el transito para conocer valores de velocidad, separation entre 
vehiculos, ocurrencia simultanea de varios vehiculos en una section, entre otros datos. 

Las conclusiones de este trabajo no son particulares para Mexico, sino que son aplicables a 
palses que no cuenten con un reglamento de pesos y dimensiones de vehiculos de autotransporte, 
un control para el cumplimiento del mismo y un modelo de carga desarrollado a partir de las 
solicitaciones reales de un pals. 

Actualmente, este grupo de trabajo sigue trabajando en esta Ilnea de investigation cuyo 
objetivo final es el desarrollo de un modelo de cargas vivas para puentes vehiculares en Mexico. 
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APENDICE A 



La clasificacion de los vehi'culos en Mexico, propuesta a principios de los setentas, consiste en 
un formato en el que, para cada pareja de caracteres, se tiene una literal que indica el tipo de 
vehiculo y un numero que senala la cantidad de ejes de dicho vehfculo. Los tipos de vehfculos y 
sus secciones son: (A) automovil, (B) autobus, (C) camion, (T) tractor, (S) semiremolque y (R) 
remolque. 



C2 

C3 

T3-S2 

T3-S3 

T3-S2-R4 



Camion de dos ejes 
Camion de tres ejes 

Tractor de tres ejes y semiremolque de dos ejes 

Tractor de tres ejes y semiremolque de tres ejes 

Tractor de tres ejes, semiremolque de dos ejes y remolque de 4 

ejes 
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RESUMEN 

El terremoto de Kocaeli, ocurrido en agosto de 1999 en Turquia, ha causado uno de los 
desastres mas graves que se han presentado en el Siglo XX. Se puede concluir de manera 
preliminar que este sismo afecto principalmente edificaciones construidas en los ultimos veinte 
anos, a lo largo de una prospera franja industrial y turistica de no mas de 15 kilometres a cada 
lado de la traza de la falla geologica de Anatolia. Miles de edificios entre cuatro y ocho pisos 
colapsaron, causando la muerte de un numero aun desconocido de personas y que podria superar 
los 30.000. Los centres urbanos mas afectados han sido Yalova, Golciik y Adapazari. En este 
ultimo, se presento la licuefaccion del suelo en un area extensa, lo que signified que muchas 
edificaciones se hundieran, se inclinaran o incluso se volcaran. Los danos en las construcciones 
se deben en su mayoria a deficiencias en la calidad de los materiales, la disposicion o detallado 
del acero de refuerzo, es decir, deficiencias en la practica de la construccion. Tambien, se pudo 
verificar que las dimensiones de las columnas y las vigas han sido insuficientes y que la mayoria 
de las configuraciones del entramado estructural no obedecen a sistemas estructurales adecuados 
para zonas sismicas, aunque Turquia ha contado desde hace varias decadas con codigos de 
construccion reconocidos a nivel internacional por su adelanto tecnico. La manera fragil como 
fallaron muchas de las estructuras ubicadas en esta zona ilustra que la existencia de disposiciones 
sismorresistentes no es una garantia de seguridad por si sola, si no se cumple y aplica dicha 
normativa como debe ser y si no hay etica en los profesionales; situacion que se repite 
aparentemente en muchos lugares del mundo. 
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ABSTRACT 

The Kocaeli earthquake, occurred recently in Turkey, produced one of the most severe natural 
desasters of the XXth century. A preliminary conclusion can be drawn that this earthquake 
mostly damaged structures built in the last twenty years along a prosperous industrial and turistic 
area extending 15 km on both sides of the Anatolian seismic fault . Thousands of buildings with 
four to eight stories collapsed, causing the death of a still unknown number of persons, which 
could exceed 30000 . The most damaged urban centers were Yalova, Golciik and Adapazari. 
In this last one, soil liquefaction occurred in an extensive area, having as a consequence 
that numerous buildings sank, tilted or overturned. But most of the buildings were 
damaged due to deficiencies in the construction practice (careless structural configuration, 
inadequate material quality, inappropriate detailing of the reinforcement), together with 
the lack of control. Although Turkey has since many decades ago building codes 
internationally recognized for their technical progress, the size of the columns and beams 
of the dammaged structures was insufficient and the greater part of the structural systems 
was not adequate for seismic areas. The fragile and sudden collapse of many of the buildings 
in the area shows that the existence of good seismic codes is not a guarantee of structural safety, 
if the code is not applied with professional ethics, fact wich occurs in many places in the world. 



1 DESCRIPCION DEL SISMO 



1.1 Localizacion geografica 

El 17 de agosto de 1999, a las 3:01 a.m. hora local (1:39.80 UTC), se presento un terremoto de 
7.4 Mw en la zona Noroeste de Anatolia (40.70° N, 20.99° E). Dicho evento causo graves danos 
en las ciudades de Yalova, Gblciik (zona epicentral), y otros pequenos poblados a Io largo de la 
costa del Mar de Marmara, pertenecientes a la provincia de Kocaeli (1.2 miliones de habitantes), 
de la cual es capital Izmit. Igualmente, fue afectada de manera muy severa la ciudad de 
Adapazari, capital de Sakarya, ubicada al Este del lago Sapanca (Figura 1). El sismo afecto 
principalmente poblaciones de una amplia zona industrial y turfstica de Turqufa, localizadas a lo 
largo de la traza de la Falla Norte de Anatolia y algunos casos de colapsos puntuales se 
presentaron en otras poblaciones mas distantes, como Estambul, aproximadamente a 90 km de la 
zona epicentral. 

El terremoto fue registrado por la red del observatorio sismologico de Kandilli y por otras 
redes a nivel mundial. La red NEIC del USGS Ie asigno a este terremoto una magnitud 7.8 Ms y 
6.3 Mb. Un numero importante de replicas se presentaron causando el colapso de un numero 
significativo de edificaciones que hablan quedado debilitadas por la accion del sismo principal. 

1.2 Sismologia y tectonica 

El epicentro del terremoto ha sido localizado all km al Sureste de Izmit, a una profundidad 
de 17 km. Esta localizacion indica que el terremoto ocurrio en uno de los segmentos de la Falla 
Norte de Anatolia; uno de los sistemas de fallamiento de rumbo lateral derecho mas estudiados y 
bien entendidos a nivel mundial en la actualidad. Se estima prehminarmente, por observaciones 
de campo de la traza de la falla en superficie, que la longitud de ruptura es aproximadamente de 
1 10 km, con un desplazamiento horizontal entre 2 y 3 m. lateral derecho, y vertical que alcanza 
en algunos sitios entre 3 y 4 m. Esta ruptura, llena una brecha o silencio sfsmico pronosticado por 
la clara migracion progresiva de sismos de Este a Oeste a lo largo del sistema de fallamiento 
desde 1939, R. S. Stein et al (1997), con once sismos con magnitudes mayores o iguales a 6.7. 
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Figura 1. Localization de epicentro y alineamiento de la Falla Norte de Anatolia 



Los terremotos sobre la Falla Norte de Anatolia son causados por el movimiento hacia el 
Norte de la Placa de Arabia contra la Placa Euroasiatica, deslizando la microplaca de Turquia 
hacia el Oeste. Desde el punto de vista sismologico, el evento tuvo una magnitud significativa 
que permite a grandes distancias lograr registros telesismicos adecuados para estudiar de manera 
confiable la directividad, la longitud de ruptura y otros parametros caracterfsticos que seguro 
seran publicados proximamente. La propagacion de replicas indica por su ubicacion que la 
ruptura de la falla va desde el Oeste de Yalova hasta el Este de Adapazari. El sismo piincipal sc 
sintio tambien en ciudades fuera de la linea de falla, por ejemplo en Ankara a 320 km. Sin 
embargo, con excepcion de algunos danos puntuales reportados en varias de estas poblaciones, 
tales como Izmit, Bursa, DUzce, Iznik, Avcilar. aparentemente debido a una atenuacion 
relativamente alta, el sismo podria decirse que solo afecto de manera significativa una franja de 
no mas de 30 km a cada lado del alineamiento de falla. 

1.3 Registros de aceleracion 

Si se tiene en cuenta que la extension de la ruptura de la falla y que la magnitud del sismo 
fueron significativas, las aceleraciones registradas durante el terremoto podrian considerarse 
sorprendemente bajas, pues no alcanzaron los 300 gal en roca o suelo duro. Sin embargo, es 
importante anotar que excepto Estambul, en las diferentes poblaciones afectadas solo se conlo 
con una estacion acelerografica. Por lo tanto no se pueden descartar aceleraciones mas severas en 
sitos especificos debido a condiciones de sitio particulares, que hubiesen podido licuar el suelo, 
como en Adapazari (407 gal), o amplificar el sismo notablemente (Figura 2). En algunas 
estaciones como las localizadas en Yarimca (322 gal) (Figura 3) y Diizce (374 gal), se evidencia 
ron efectos de campo cercano y directividad que caracterizan un evento como el descrito. 











26 



O. Cardona y A. Barbat 




Figura 2. Valores de aceleracion en %g registrados en la zona afectada. 



Es importante considerar que el contenido frecuencial del sismo pudo haber tenido una 
influencia importante en la respuesta de las estructuras cuyos periodos de vibracion se 
encontraban entre 0.5 y 1.0 seg, lo que explicarfa en parte la concentracion de danos y colapsos 
de edificaciones de cuatro a ocho pisos, que son la mayorfa en la zona afectada (Figuras 3 y 4). 
Esto tambien explicarfa porque edificaciones rfgidas de uno y dos pisos, de periodos de vibracion 
menor, no se afectaron de la misma forma. Este hecho sorprende, dado que un numero importante 
de edificaciones de mamposterfa no reforzada de ladrillo o de adobe de tierra con refuerzos de 
madera no sufrieron dafio. 



STATION YPT ■ NS 




STATION YPT - EW 




STATION YPT - UD 




Figura 3. Registros de aceleracion horizontal y vertical en la Estacion Yarimca. 
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Figura 4. Espectro de respuesta para el registro de aceleracion en Yarimca. 

Per tratarse de un sismo de campo cercano para edificaciones ubicadas cerca de la zona 
epicentral, la directividad y la energfa sfsmica aparentemente contribuyeron de forma severa a 
que las edificaciones de cuatro o mas pisos fallaran de manera fragil a nivel del primer piso. Una 
muy alta demanda de energfa inelastica, que no fue posible atender por parte de las estrucuturas, 
explicarfa porque muchas edificaciones colapsaron perdiendo su primera planta y en algunos 
casos la segunda, sin presentar mayores danos en los pisos superiores. 

1.4 Efectos de sitio 

El caso mas evidente de efectos locales de caracter geotecnico, fue la licuefaccion de suelos. 
En Adapazari una extensa area de la zona urbana se licuo por la accion del sismo, causando el 
hundimiento de no menos de 50 cm en promedio de cientos de edificios (Foto 1). Muchos de 
el los se inclinaron debido a asentamientos diferenciales e incluso hubo un amplio numero de 
casos de vuelco de edificaciones cuya cimentacion era una losa relativamente superficial (Foto 
2). Las tuberfas de agua se rompieron debido a que las calles no se hundieron igualmente y por 
esta razon los primeros pisos se inundaron con aguas de las redes. Este fenomeno se presento en 
forma generalizada en Adapazari debido a que la ciudad se encuentra sobre suelos de depositos 
aluviales arenosos sueltos y saturados hasta profundidades de 5 m, donde el nivel freatico se 
encuentra muy superficial, entre 1 y 3 m, por la cercanfa de dos rfos y el Lago Sapanca. Estos 
suelos, en ocasiones se encontraban sobre otros suelos de grava y arena mas densos o debajo de 
llenos artificiales de poco espesor, A. Erken (1999). 

Tambien se reportaron casos de licuefaccion en Izmit y Yalova, cerca a la tinea de costa, 
donde se presentaron hundimientos, inclinaciones y colapsos de edificaciones. En estas 
localidades, los suelos arenosos sueltos se encontraban subyaciendo depositos delgados de 
arcilla, e igualmente la tabla de agua se encontraba aproximadamente a 1 m de profundidad. 
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Foto 2. Edificio volcado por licuefaccion del suelo y aceleracion intensa en Adapazar. 
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Otros efectos de sitio, tales como amplificacion por suelo blando y amplification topografica, 
parecen evidentes en lugares como Avcilar, cerca a Estambul, y Diizce, donde las aceleraciones 
registradas asi lo indican. Sin embargo, el desconocimiento de los tipos de suelos en la mayor 
parte de la zona afectada y la falta de registros acelerograficos no perrwten deducir con certeza 
que tan generalizada pudo ser la amplificacion sismica. Queda la duda, no obstante, si muchos 
casos de colapso, un tanto incomprensibles, en muchos casos se pudieran explicar por 
condiciones de sitio desfavorables que pudieron amplificar notablemente el movimiento. 

La expresion en superficie de la ruptura de la falla causo un numero notable de colapsos y 
danos debido a las deformaciones excesivas el suelo tanto en el sentido horizontal como vertical 
(Foto 3). En su mayoria, los edificios que fueron cruzados por la traza o que se encontraban muy 
cerca a ella se afectaron gravemente o colapsaron (Fotos 4 y 5). 




Foto 3. Traza en superficie de la falla de rumbo lateral derecho, USGS. 



1.5 Sismicidad y amenaza sismica 

Esta zona ha sido afectada frecuentemente por terremotos que han causado danos. Existen 
registros historicos bien documentados de diez y ocho sismos que desde el ano 29 DC han 
afectado de manera significativa la zona, de los cuales hay descripciones detalladas de sus 
efectos. Durante ei presente siglo, se recuerdan los sismos del 26 de mayo de 1957 y 22 de julio 
de 1967 que afectaron especi'ficamente la misma zona afectada por el terremoto del 17 de agosto 
de 1999. Otros eventos que han ocurrido en las vecmdades son el del 20 de junio de 1943 (Ms 
6.4) al Este de Adapazari, en Hendek, y del 18 de septiembre de 1963 (Ms 6.4) en Valova 





Foto 4. Polideportivo en Golciik, atravesado por la traza de la falla. 



Foto 5. Cimentacion del polideportivo elevada por subsidencia del suelo. 
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Desde los anos 70 se han realizado para Turqufa un amplio numero de estudios de amenaza o 
peligrosidad sfsmica, que hace esta una de las zonas del mundo mas estudiadas. Considerando 
catalogos sfsmicos completos, antes de 1963 para M > 5.5 y despues de 1963 para M > 4, los 
estudios mas recientes, M.Erdik, et al (1999), que tienen en cuenta leyes de atenuacion cada vez 
mas precisas, indican que la zona donde ocurrio el sismo tiene aceleraciones pico del suelo del 
orden de 0.5g con un 10% de excedencia en 50 anos (Figura 5). Para efectos de diseno 
sismorresistente, esta zona se encuentra en la Zona 1 (0.4g) de cuatro establecidas, en la cual las 
exigencias del codigo de construcciones son las mas rigurosas, Ministry of Public Works and 
Settlement of Turkey (1999). 



EARTHQUAKE ZONING MAP FOR TURKEY j 




Provincial * 
Capitals 



Figura 5. Mapa de zonificacion sfsmica de Turqufa para aplicacion del codigo. 



2 EFECTOS DEL TERREMOTO 
2.1 Dafios en edificaciones 

Los dafios se concentraron en su mayorfa en edificaciones relativamente modemas, 
construidas en los ultimos veinte anos. El desarrollo de la region en ese perfodo ha sido muy 
notable debido a la polftica de industrializacion y estfmulo al turismo, lo que se ha traducido en 
un crecimiento muy importante de la poblacion y la construccion acelerada de edificaciones que 
se caracterizan en su mayorfa por ser estructuras aporticadas de concreto reforzado de 4 y 5 pisos 
y en algunos casos hasta de 8 pisos, a lo largo de la costa del Mar de Marmara y en las cercanfas 
del Lago Sapanca (Fotos 6 y 7). 

Curiosamente, a diferencia de otras zonas del mundo, no existe en el area afectada por el 
terremoto un uso extensivo de edificaciones de mamposterfa no reforzada de ladrillo o adobe, 
dado que el poblamiento de la zona es muy reciente. Las muy pocas edificaciones de 
mamposterfa no reforzada que se encontraron, sorprendentemente, no presentaron danos o fueron 
menores. A continuacion se examinan las principales causas de los danos en los edificios. 
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Foto 6. Vista de edificios de varios pisos colapsados en la zona afectada. 
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2.1.1 C alidad de los materiales 

Un amplio numero de edificios se construyeron con materiales de muy mala calidad. El acero 
de refuerzo utilizado practicamente en la totalidad de las estructuras era acero liso no ductil de 
baja resistencia (2,100 kg/cm 2 ), con una capacidad de elongacion que posiblemente no alcanza a 
superar el 5%; valor muy por debajo del exigido como mlnimo para el acero en zonas slsmicas, 
por las normas actualmente en diferentes partes del mundo (Foto 8). Este aspecto de especial 
importancia para garantizar el comportamiento ductil y de alta capacidad de disipacion de energla 
inelastica, fue el gran ausente en este terremoto, posiblemente debido a que la introduccion del 
acero corrugado de alto porcentaje de elongacion ha sido reciente. La falla fragil del acero fue 
posible identificarla en un amplio numero de casos donde se podia ver sin dificultad el refuerzo 
expuesto de los elementos estructurales. Por otra parte, la corrosion del acero ha sido 
generalizada en la zona por la falta de recubrimiento debido a las dimensiones insuficientes de 
los mismos. Se pudo detectar no solo en los edificios danados sino en estructuras en obra negra 
en proceso de construccion, que el acero en ocasiones queda expuesto, sin recubrimiento de 
hormigon alguno, en una zona que ambientalmente por su contenido de salinidad no podrla 
esperarse otra cosa que un accion corrosiva acelerada (Foto 9). Tambien, es importante indicar 
que solo se observaron varillas de diametros menores, lo que es tlpico con aceros lisos y de baja 
resistencia (Foto 10). 

A primera vista tambien fue posible detectar que la calidad del hormigon es bastante pobre. 
Solo a principios de los anos 80 se inicio el uso de hormigon premezclado el cual se usa actual- 




Foto 8. Nudo de viga-columna con acero liso insuficiente y de baja ductilidad. 
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mente en un 40% de las construcciones. Es decir, en la mayoria de los casos el hormigon se 
realiza en obra y su control de resistencia es muy deficiente. No solo los valores de resis-tencia 
son bajos sino que su dispersion es muy alta. Infortunadamente, en muchas edificacio-nes en 
proceso de construccion se pudo comprobar que los elementos estructurales quedan mal vibrados 
y que el hormigon no alcanzaba a Uegar a todos los sitios del encofrado; que-dando asi con zonas 
de mayor concentracion de agregado y poco cemento. Esto posiblemente se debe a que las 
dimensiones de las columnas y las vigas son muy insuficientes y a que no se utilizan tecnicas 
adecuadas para la vibracion del hormigon una vez colocado (Foto 1 1 y 12). 




Foto 11. Estructura en construccion donde se detecta la falta de recubrimiento. 



2. 1.2 Control y supervision 

Se puede deducir que las edificaciones afectadas, en su mayoria. se han disenado y construido sin 
un control de calidad adecuado y un debido seguimiento o supervision de la obra. Existe un 
consenso en los especialistas locales que el control no ha sido el comun denominador y que 
incluso es preocupante que las edificaciones se construyan con mano de obra muy deficiente en 
terminos tecnicos. Aparentemente, la posibilidad de poder cubrir con acabados arquitectonicos 
las estructuras que claramente dejan ver en obra negra que su construccion ha sido muy 
deficiente, se ha convertido en la altemativa poco etica de algunos profesionales para evitar 
demoliciones o reparaciones de elementos estructurales que en la construccion no han quedado 
adecuadamente. De acuerdo con las indagaciones realizadas, no se realiza ningun tipo de prueba 
de la resistencia del hormigon ni del acero ni se hacen otras pruebas clasicas o ensayos que se 
exigen en las normas o reglamentos de construccion en relacion con la supervision tecnica. 
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Foto 12. Estructura en construction donde se aprecia la falta de vibrado del hormigon. 



2.1.3 Configuration estructural 

Las edificaciones afectadas por el terremoto en su mayoria no cumplian con el codigo de 
construction sismorresistente. Turqufa ha tenido un codigo de construcciones sismorresistentes 
desde 1940, que se ha actualizado tecnicamente a lo largo de los anos. La ultima actualization se 
realizo oficialmente en julio de 1998 por el Ministerio de Obras Publicas y Asentamientos 
Humanos, por lo tanto, no ha existido un descuido en cuanto mantenerse al dia en materia de 
sismorresistencia. Sin embargo, es totalmente evidente que dicha normativa no se ha aplicado 
debidamente. Las disposiciones actuales como posiblemente las anteriores indican. por ejemplo. 
que la dimension minima de una columna no puede ser menos de 25 cm y que el area minima de 
la section no debe ser inferior a 750 cm 2 , situation que no se cumple en ningun caso en las 
edificaciones construidas en la zona. Es comun encontrar columnas en los cuales el ancho de la 
section rectangular tiene 15 cm o en el mejor de los casos 20 cm, posiblemente con el fin de que 
las columnas queden de un espesor similar al de las paredes o tabiques divisorios (Foto 13 1 . 

Esta circunstancia indica que las columnas son excesivamente esbeltas, mas si es comun 
encontrar que la primera planta de los edificios frecuentemente es de mayor altura que las demas. 
posiblemente con el fin de promover su uso para comercio. Ahora bien, el otro lado o cara de la 
section en las columnas usualmente es de 40 cm o 50 cm, lo que da una apariencia de notable 
rigidez, particularmente en edificaciones de dos y tres pisos, sin embargo la debilidad y esbeltez 
de la columna en la otra direction ha contribuido notablemente a que muchas edificaciones 
fallaran por insuficiencia de rigidez. Este hecho lo agravo, en muchos casos, que los entrepisos 
fueran losas planas o placas aligeradas muy delgadas (Foto 13), lo que signified que 
contribuyeron a dar una mayor flexibilidad y por lo tanto mayor deformation o deriva ante careas 
laterales, dado que no hubo una action de portico resistente a momento. En consecuencia se 
puede firmar que hubo innumerables casos de inestabilidad local o global por efecto P-Delta. 
Tambien ha sido generalizado el uso de losas de concreto reforzado vaciadas en el sitio. de 
maximo 10 cm entre vigas y 5 cm en voladizos de balcones, con vigas en algunos casos 
relativamente alias de 40 cm, propensas al pandeo lateral debido a que el ancho de la section no 
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superaba en ningun caso los 15 cm. El pandeo lateral se registro debido a la debilidad de la union 
en los nudos con las columnas a las cuales las vigas llegaban en forma excentrica frecuentemente 
(Fotos 14 y 15). Por otra parte, es importante destacar que en la mayorfa de los casos, las 
columnas no conformaban con las vigas porticos bien definidos en un piano, es decir que las 
vigas en ocasiones se apoyan sobre otras vigas y que las columnas no se encuentran alineadas. 
Aparentemente, estos aspectos parecen ser una practica comun que se ha replicado en la 
construccion informal a pesar de que el codigo de construcciones desestimula este tipo de 
practicas, castigandolas como irregularidades de configuracion para el analisis y diseno. Es 
importante mencionar que un numero significativo de edificaciones les fallo el primer piso por la 
ausencia de muros que evitaran la accion del efecto de piso debil o blando y en muchos casos por 
la perdida de muros usualmente construidos sin conexion alguna a la estructura y con ladrillo de 
perforation horizontal que es altamente fragil o quebradizo. 




Foto 13. Elementos estructurales de dimensiones insuficientes. 



Esta circunstancia indica que las columnas son excesivamente esbeltas, mas si es comun 
encontrar que la primera planta de los edificios frecuentemente es de mayor altura que las demas, 
posiblemente con el fin de promover su uso para comercio. Ahora bien, el otro lado o cara de la 
seccion en las columnas usualmente es de 40 cm o 50 cm, lo que da una apariencia de notable 
rigidez, particularmente en edificaciones de dos y tres pisos, sin embargo la debilidad y esbeltez 
de la columna en la otra direction ha contribuido notablemente a que muchas edificaciones 
fallaran por insuficiencia de rigidez. Este hecho lo agravo, en muchos casos, que los entrepisos 
fueran losas planas o placas aligeradas muy delgadas (Foto 13), lo que signified que 
contribuyeron a dar una mayor flexibilidad y por lo tanto mayor deformation o deriva ante cargas 
laterales. dado que no hubo una accion de portico resistente a momento. En consecuencia se 
puede firmar que hubo innumerables casos de inestabilidad local o global por efecto P-Delta. 





38 



O. Cardona y A. Barbat 




Foto 15. Mala configuracion estructural de entramados que forman porticos 
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Tambien ha sido generalizado el uso de losas de concrete reforzado vaciadas en el sitio, de 
maximo 10 cm entre vigas y 5 cm en voladizos de balcones, con vigas en algunos casos 
relativamente altas de 40 cm, propensas al pandeo lateral debido a que el ancho de la section no 
superaba en ningun caso los 15 cm. El pandeo lateral se registro debido a la debilidad de la union 
en los nudos con las columnas a las cuales las vigas llegaban en forma excentrica frecuentemente 
(Fotos 14 y 15). Por otra parte, es importante destacar que en la mayorla de los casos, las 
columnas no conformaban con las vigas porticos bien definidos en un piano, es dec ir que las 
vigas en ocasiones se apoyan sobre otras vigas y que las columnas no se encuentran alineadas. 
Aparentemente, estos aspectos parecen ser una practica comun que se ha replicado en la 
construction informal a pesar de que el codigo de construcciones desestimula este tipo de 
practicas, castigandolas como irregularidades de configuration para el analisis y diseno. Es 
importante mencionar que un mimero significativo de edificaciones les fallo el primer piso por la 
ausencia de muros que evitaran la action del efecto de piso debil o blando y en muchos casos por 
la perdida de muros usualmente construidos sin conexion alguna a la estructura y con ladrillo de 
perforation horizontal que es altamente fragil o quebradizo. 

Aunque se podria decir que algunos edificios podrfan ser porticos con tabiques de ladrillo, o 
incluso de mamposterla de ladrillo confinada, debido a la disposition arquitectonica de las 
columnas, en muchos casos, los muros en la periferia del edificio no quedan confinados por los 
porticos, lo que significa que las paredes o tabiques no solamente pueden caerse facilmente sino que 
no contribuyen aportar rigidez a la estructura (Foto 16). Una practica muy perjudicial, 
aparentemente promovida por las autoridades pollticas y sin ningun control, ha sido el construir 
formal o informalmente nuevos pisos superiores sobre estructuras existentes, lo que cambia 
obviamente el comportamiento de la estructura original que desde su inicio, podria decirse, es 
deficiente para atender las solicitaciones de las cargas laterales causadas por sismos. Tambien, en el 
caso de estructuras prefabricadas de concrete reforzado, usualmente construidas para naves 
industrials, se pudo detectar graves deficiencias en las conexiones. La falta de un diafragma a nivel 
de la cubierta facilito que estas edificaciones cuando no contaban con arriostramientos laterales, 
tales como muros, colapsaran facilmente. 




Foto 16. Mamposterla de ladrillo hueco sin conexion adecuada en la periferia. 
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2.1.4 Disipacion de energia 

El desconocimiento o no la aplicacion del codigo de construcciones es aun mas evidente en la 
falta de un adecuado detallado del refuerzo. La insuficiencia del acero y su colocacion 
descuidada en las columnas, vigas y nudos de las estructuras son sin duda la causa de la falla 
fragil y subito generalizado en las edificaciones afectadas. La mayoria de los colapsos parciales o 
totales se debieron a un pesimo comportamiento de las estructuras que no tuvieron un 
comportamiento ductil que les permitiera disipar energia inelastica. La falta de acero transversal 
en vigas y columnas ha sido la constante en este tipo de edificaciones, donde por falta de 
confinamiento cerca de los nudos no fue posible que se desarrollara una capacidad en terminos de 
resistencia y disipacion de energia ante las solicitaciones sismicas. Los flejes o estribos eran 
insuficientes y usualmente se encontraban muy separados (Fotos 17 y 18). El gancho no doblaba 
a 135° sino a 90° lo que facilito que los aros se abrieran y el acero longitudinal se pandeara. De 
igual forma, se puede afirmar que el acero longitudinal en su mayoria fue insuficiente si se tiene 
en cuenta que en muchas edificaciones es acero liso de baja resistencia. Errores de despiezo se 
encontraron en muchas estructuras, en las cuales no se logro un adecuado desarrollo del refuerzo 
por falta de anclaje y adherencia (Foto 19). 




Foto 17. Columna con poco confinamiento de acero transversal. 



Uno de los aspectos de mayor incidencia eh la falla fragil y subita de las edificaciones ha sido 
la colocacion precaria del acero en los nudos y la deficiencia de los mismos por su debilided 
inherente, debido a la excentricidad de las vigas en su conexion con las columnas y por la falta de 
dimensiones adecuadas de los elementos estructurales (Fotos 20, 21, 22 y 23). En general, puede 
afirmarse que no hubo la capacidad de deformacion sin falla y que las estructuras failaron 
subitamente sin disipar energia inelastica. Por esta razon, se presento una situacion poco mual 
despues de la ocurrencia de un sismo: las edificaciones en su mayoria presentaron dahos seseros. 
que se traducen en muchos casos en un colapso, o no presentaron danos estructurales. Solo en 
muy pocos casos se encontraron danos graves, moderados o leves, tipicos de una respuesta 
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Foto 21. Degradacion de un nudo de viga-columna y articulacion de la estructura. 
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Foto 23. Edificio articulado en primera planta por mecamsmo de rotulas plasticas. 



Es importante mencionar, no obstante, que en lugares cercanos a la zona epicentral donde el 
remoto en terminos de ingenieri'a podrfa considerarse de campo cercano, la directividad y el 
■acter imnulsivo el sismo tambien'contribuyeron a que no se alcanzara a desarrollar un 
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mecanismo de disipacion de energfa, incluso en edificaciones competentes. El diagrama de 
radiacion del sismo indica que las dos componentes de aceleracion horizontal registradas 
coinciden en ser las mas fuertes aproximadamente en el mismo instante, lo que significa un fuerte 
impacto sobre las edificaciones que impide desarrollar ciclos de histeresis resultado de la accion 
altemada de esfuerzos debido a la vibracion u oscilacion de la estructura (Figura 6). 




Figura 6. Ilustracion de la simultaneidad de las componentes sismicas horizontales. 
EstacionYPT, terremotode Kocaeli, 17.08.1999. 



2. 1. 5 Informalidad de la construccion 

A pesar de contar con una norma sismorresistente idonea, su desconocimiento, la falta de 
control y en algunos de etica profesional han contribuido significativamente a que las 
edificaciones sean altamente vulnerables. Tal como se menciono antes, por emulacion en la zona 
se han construido un numero muy amplio de edificaciones sin la participacion de ingenieros y 
arquitectos conocedores del codigo y de los requisitos en general de comportamiento 
sismorresistente. Estas edificaciones que a primera vista se les ve con una adecuada rigidez. 
usualmente solo la tienen en una direccion, siendo muy insuficiente en la otra, dado que las 
columnas son muy delgadas y esbeltas. Sin duda la practica de la construccion informal de 
edificaciones de concreto reforzado no ha sido buena. Lamentablemente, la informalidad en la 
construccion en esta zona se ha fomentado incluso por politicas que han promovido la 
construccion de varios pisos adicionales sobre las edificaciones existentes. Este aspecto 
seguramente contribuyo al mal comportamiento de muchas de ellas. 

No deja de ser asombroso, sin embargo, que edificaciones antiguas de mamposteria no 
reforzada de ladrillo, adobe y madera, incluso mesquitas y sus esbeltos minaretes no presentaran 
danos severos. Algunas de ellas se le ve en muy buen estado. Adicionalmente, es posiblemente 
que el contenido frecuencial del sismo las haya favorecido. Algunos espectros de respuesta 
indican que la aceleracion en edificaciones con perfodos entre 0.4 y 0.7s pudo llegar a ser del 
orden del 0.7g, lo cual podria explicar la concentration de danos en el rango de las edificaciones 
de concreto reforzado afectadas por el terremoto. 

Una de las preocupaciones mayores a nivel tecnico es la falta de una certification de 
idoneidad profesional para el ejercicio en el campo del analisis y diseno estructural > la 
geotecnia. Dada la proliferation de facultades de ingenierfa algunos creen que el nivel de los 
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profesionales ha bajado notablemente en los ultimos anos y que la experiencia de los disenadores 
y constructores no es suficiente para emprender proyectos con la debida seguridad y 
confiabilidad. Seguramente este problema tendra que ser considerado en conjunto con la 
necesidad de lograr un mayor control y vigilancia no solamente de la idoneidad de quien disena o 
construye, sino tambien si quien construye es un profesional o no. 



2.1.6 Planificacion urbana 

La no consideracion de los peligros naturales en el ordenamiento territorial y la planeacion de 
los usos del suelo urbano ha sido tambien un aspecto que ha contribuido notablemente a 
extension de los danos. La construccion sobre suelos licuables y sobre trazas de la falla en 
superficie es una circunstancia que se podrfa evitar o desestimular, sin embargo no ha sido 
debidamente considerada como un determinante del desarrollo urbano. En este desastre fue 
evidente que no es posible evitar los danos en las edificaciones cuando estas son afectadas por 
deformaciones excesivas del suelo. La licuefaccion del suelo afecto areas urbanas extensas y la 
expresion en superficie de la falla causo el colapso de innumerables edificaciones. Algunos 
consideran que el proceso de descentralizacion administrativa ha causado que la planificacion 
regional que antes hacia el gobiemo central con recursos humanos idoneos, se realice 
deficientemente o sencillamente no se realice. Dado que crecimiento de la poblacion y la 
construccion de edificaciones ha sido muy acelerada en los ultimos anos en la zona, se podrfa 
deducir que incluso para el gobiemo central hubiese sido muy diffcil controlar el asentamiento 
humano en lugares inadecuados, mas si se sabe que existen fuerzas economicas y sociales 
relacionadas con la tenencia de la tierra que no son faciles de neutralizar. Se espera que esta 
experiencia sirva para que se reoriente el desarrollo urbano y se tengan en cuenta los efectos 
sfsmicos colaterales antes mencionados en la planificacion de los futuros programas de 
reconstruccion. 

2.2 Danos en la industria 

La zona afectada por el terremoto se caracteriza por su rapido desarrollo industrial en los 
ultimos anos. Allf se encuentran plantas petroqufmicas, ensambladoras de automoviles, fabricas 
de papel, plantas de produccion de acero, cemento, etc. La industria petroqufmica de la zona del 
Mar de Marmara fue la mas afectada por el terremoto. Se presentaron incendios en la refinerfa 
Tupras en Korfez que duraron vario dfas, en las afueras de Izmit. Danos importantes se 
registraron en instalaciones, torres, gruas y tanques que fallaron por a accion sfsmica. Los costos 
de los danos en el sector se estiman en US$ 2,000 millones (2 billones) de dolares, que en su 
mayor parte se encontraban asegurados. Los danos en la industria, tanto estructurales como no- 
estructurales seguramente seran objeto de estudios especiales debido a las implicaciones de lucro 
cesante que significan, particularmente debido a que las aceleraciones no fueron muy fuertes. 

2.3 Danos en lineas vitales 

Los servicios publicos en la mayorfa de las poblaciones del area fueron gravemente afectados. 
Un mes y medio despues del sismo aun permanecen algunos sitios sin agua potable, sin energfa y 
telefonos en la zona por danos locales. La central de 380 kV de Adapazari quedo fuera de 
servicio pero pocos dfas despues del evento se puso en funcionamiento de nuevo. Debido a la 
licuefaccion del suelo las tuberfas de agua se rompieron y cambiaron de pendiente en sitios como 
Adapazari. En relacion con las vfas, el paso estuvo restringido solo por pocos dfas mientras se 
llevaron a cabo las reparaciones respectivas. El movimiento de la falla causo la interrupcion del 
trafico de vehfculos en varfas vfas, particularmente en Arifiye en donde colapsaron dos puentes 
por perdida de apoyo del tablero sobre los estribos y la traza de la falla rompio la calzada en 
varios sitios (Foto 24). 
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Foto 24. Traza de la falla deformando el suelo excesivamente de una via vehicular. 



Igualmente, el trafico ferroviario se afecto notablemente por el daiio de los rieles. Varios 
puertos sobre el Golfo de Izrnit se afectaron por licuefaccion o por la deformacion excesiva del 
suelo causado por la falla; el mas importante de los puertos afectados fue la Base Naval en 
Gblciik. Algunos lugares de la costa bajaron varios metros razon por la cual algunas calles, 
parques y muelles de embarcaciones menores se hundieron y fueron inundados por el mar. 

2.4 Impacto social y economico 

De acuerdo con la informacion oficial del General Directorate of Disaster Affairs Earthquake 
Research Department of Turkey, las cifras de perdidas indican que 15,135 personas perdieron la 
vida, pero se desconoce el numero de desaparecidos que se estima puede ser una cifra del mismo 
orden de magnitud. Se atendieron 24,000 heridos en los hospitales. El numero de afectados se 
calcula en 600,000, que perdieron su vivienda y sus espacios de trabajo. El terremoto causo danos 
moderados, graves y perdida total a cerca del 20,000 unidades de vivienda que no son reparables, 
de las cuales 50,000 sufrieron danos graves o colapsaron. Se estima que 2,000 edificios 
colapsaron y 4,000 sufrieron danos severos lo que puede estimarse que significa una perdida 
economica directa de aproximadamente US$ 5,000 millones de dolares (US$ 5 billones). Los 
danos en la industria se estiman del orden de US$ 2,000 millones de dolares. Los danos en vlas y 
puertos se calculan en US$ 1,400 millones de dolares y en general se estima que incluyendo otras 
perdidas de caracter socieconomico, las pedidas pueden sumar US$ 16,000 millones de dolares. 
lo que significa un impacto total del sismo del orden del 7% del PIB del pals. Este terremoto 
afecto gravemente una de las zonas mas prosperas de Turquia, de notable desarrollo industrial en 
los ultimos veinte anos y con un importante desarrollo turfstico que se tradujo en una notable 
actividad de la construccion de edificaciones en los ultimos anos. Aun cuando los danos en la 
industria se encuentran aparentemente cubiertos en su totalidad, de las edificaciones afectadas. se 
estima que solo 8% estaban cubiertas por seguros. 

Este evento compromete de manera severa el desarrollo economico y social de 
Turquia, que en los ultimos anos ha realizado esfuerzos para lograr pasar a la categona de 
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pafses desarrollados. Infortunadamente, el impacto en la industria y el proceso de 
reconstruccion implicaran probablemente gastos que cambiaran las polfticas y programas 
de gobierno y su nivel de endeudamiento, lo que se traduce en empobrecimiento, retrazo 
por varios anos y un obstaculo mas para lograr entrar en el primer mundo. 



3 ENSENANZAS 

Del desastre de Kocaeli en Turquia, quedan varias ensenanzas importantes para la ingenierfa 
si'smica y la prevencion de desastres: 

• Se demostro de nuevo el peligro que representa para una region o un pafs el contar con 
un inventario de edificaciones que no cumplen con requisitos mfnimos de 
simorresistencia. Las edificaciones afectadas por este terremoto en su mayoria han sido 
estructuras de concrete reforzado con elementos estructurales de dimensiones deficientes, 
mala calidad del hormigon e insuficiente acero de refuerzo. 

• Los principales danos en las estructuras se presentaron por esfuerzo cortante, debido al 
deficiente detallado del acero de . refuerzo y la falta de acero transversal y de 
confinamiento del hormigon en los elementos estructurales. Las fallas fueron de tipo 
fragil, lo que demostro una deficiente capacidad de disipacion a absorcion de energfa 
inelastica. 

• Una leccion muy importante para otros lugares del mundo, es que es necesario mas no es 
suficiente con tener un codigo de construcciones sismorresistentes actualizado de 
acuerdo con el estado del conocimiento. Debe existir ademas un control y supervision 
que garantice que la normativa se aplica correctamente. 

• Se ratified la necesidad de tener en cuenta los fenomenos sfsmicos colaterales, como la 
licuefaccion del suelo y las zonas de ruptura de fallas geologicas, en los planes de 
ordenamiento urbano y reglamentacion de usos del suelo, con el fin de evitar que estas 
zonas sean urbanizadas. 

• Una vez mas, se ratified la importancia de contar, a nivel de las ciudades, con programas 
i nteri nstitucionales de mitigacion de riesgos que involucren aspectos tecnicos y 
cientificos, de planificacion y desarrollo, de educacion e informacion publica y de 
preparacion para atender emergencias y desastres. 

Muchas regiones y pafses del mundo tienen situaciones similares. Edificaciones con las 
mismas caracterfsticas que las hacen altamente vulnerables se han construido y se siguen 
construyendo en muchas ciudades del mundo. En esta ocasion le ha tocado a Turquia, pero 
probablemente en los proximos anos esta situacion se repetira en diversos lugares. 
Infortunadamente, para que haya un desastre sfsmico no solo se requiere de que exista una alta 
amenaza o peligrosidad sfsmica, sino una alta vulnerabilidad de las edificaciones, lo cual 
usualmente esta asociado con el hecho de que la amenaza sfsmica se considere poco importante. 



4 AGRADECIMIENTOS 



Este informe se realizo con base en la experiencia directa de visitar la zona de desastre entre el 
20 y el 27 de octubre de 1999, visita financiada por la empresa Calidad Siderurgica de Espafia.. 




48 



O. Cardona y A. Barbat 



Parte de la informacion aquf presentada ha sido facilitada por Istanbul Technical University, Civil 
Engineering Faculty, ( http://Mantle.geopdtu.edu.tr/') Istanbul; Bogazici University, Kandilli 
Observatory ( http://www.koeri.boun.edu.tr/earthQk/earthquake.htTn ). Istanbul; General 
Directorate of Disaster Affairs, Earthquake Research Department, ( http://www.deprem.gov.tr /): 
Turkey, Earthquake Engineering Research Institute (http://www.eeri.org '): United States 

Geological Survey (http://groundmotion.cr.usgs.gov/izm/) : y The Southern California 
Earthquake Center (SCEC) (http://www.scec.org/research/turkev.html ). 

Especial reconocimiento se desea expresar a los profesores Esin Inan, Faruk Karadogan, Atilla 
M. Ansal, Ayfer Erken de Istanbul Technical University, a los profesores Mustafa Erdik y Ozal 
Yuzugullu de Bogazici University Kandilli Observatory and Earthquake Research Institute. 
Igualmente, se agradece al Dr. Konuralp Girgin y la Ing. Pelin Gundes Bakir por haber 
colaborado en la visita a las aeas afectadas. 



5 REFERENCIAS 



R. S. Stein, A. A. Barka y J. El. Dieterich, 1997, “Progressive failure on the North Anatolian fault 
since 1939 by earthquake stress triggering”. Geophysical Journal International, 1997 , 128, 594- 
604. 

http://quake.wr.usgs.gOv/study/turkev/#background 

M. Erdik, Y. Alpay Biro, T. Onur, K. Sesetyan y G. Birgoren, 1999, “Assessment of earthquake 
hazard in Turkey and neighboring regions”, Bogazici University Kandilli Observatory and 
Earthquake Research Institute, Istanbul. 

http://www.koeri.boun.edu.tr/earthqk/earthauake.htm 

A. Erken, 1999, “The effect of soil condition during Kocaeli earthquake”, Intanbul Technical 
University, Civil Engineering Faculty, Maslak, Istanbul. 
http://Mantle.geop.itu.edu.tr/ . 

Ministry of Public Works and Settlement, 1999, “Specification for Structures to be Built in Disaster 
Areas: Part III - Earthquake Disaster Prevention”, (Chapter 5 through Chapter 13), Issue 
2.9.1997, Official Gazette No.23098, effective from: 1.1.1998, amended on: 2.7.1998, Official 
Gazette No.23390, Government of Republic of Turkey. 




Revista internacional de 



Vol. 4-2, 5-1, 49-60 (2000) 



Ingenieria 
de estructuras 



Modelacion de la atenuacion espacial de 
la energia sismica por medio de redes 
neuronales 



Jorge E. Hurtado y Manuel Builes 

Universidad Nacional de Colombia 
Apartado 127, Manizales, Colombia 



Resumen 

La estimation de la atenu'acion espacial de la energia sismica constituye un paso impor- 
tante en los estudios de vulnerabilidad y dario causados por sismo dado en lugares donde 
no se hayan obtenido registros. En este articulo se explora la posibilidad de emplear los 
dispositivos de inteligencia artificial conocidos como redes neuronales multi- estrato para el 
calculo de los contornos de iso-aceleracion a partir de los datos instrumentales disponibles 
en algunos puntos debido a la habilidad de tales redes para modelar sistemas altamcnte no 
lineales y complejos. El estudio se realiza sobre los datos de un sismo ocurrido en Colombia 
y registrado por la red sismica nacional. Los resultados muestran una correlation favorable 
con las curvas de igual intensidad de Mercalli observada, por lo que se puede concluir que 
tales dispositivos constituyen una tecnica con amplias perspectivas para dicho proposito. 



Summary 

The assessment of the spatial attenuation of the seismic energy is an important task 
to perform in calculating the vulnerability and damage caused by a given earthquake in 
those places where there are no records available. In this paper, the possibility of using 
an artificial intelligence device, known as multi-layer neural networks for drawing equal 
acceleration contours on the basis of a set of instrumental data points is explored due to their 
ability for modelling highly nonlinear and complex systems. The study is performed upon 
the information of an earthquake occurred in Colombia recorded by the Colombian seismic 
network. It is shown that the estimations afforded by the neural networks compare well 
with the isoseismals map of the event, so they are a promising technique to the mentioned 
purpose. 
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1. INTRODUCTION 

En la actualidad existen muchas medidas de la intensidad sismica y de sus efectos desde el 
punto de vista de la ingenieria. Entre ellas pueden citarse los valores pico de la aceleracion, 
la velocidad y el desplzamiento, la aceleracion espectral y la intensidad de Mercalli asi 
como diversas medidas de la energia sismica. Todas estas variables presentan un patron 
de radiacion desde el epicentro en el cual se observa que ellas disminuyen con la distancia. 
Sin embargo, la forma de esta atenuacion varia de evento a evento. Esto se debe a una 
combinacion de factores tales como la orientacion de la zona de ruptura de la falla, la 
geologia regional, la topograffa, etc. Como consecuencia de esto, el estudio de la atenuacion 
espacial debe realizarse para cada evento, sin que sea viable la utilization de los resultados 
de estudios anteriores similares. 

Los contornos de iso-aceleracion se encuentran mucho menos estudiados que los de 
intensidad de Mercalli. En principio, se sabe que la forma de estos ultimos depende 
basicamente de la transmisibilidad de las ondas en las dos direcciones, la cual suele presentar 
una marcada anisotropia. Por otra parte, los contornos tienden a seguir las tendencias 
orograficas y tectonicas (Newmark y Rosenblueth, 1976). La importancia de estudiar los 
contornos de iso-aceleracion resulta patente si se considera que, en multiples casos, se 
requiere estimar los va’lores indicatives mencionados anteriormente en sitios en los cuales 
no se ha obtenido ningun registro instrumental. Tal es el caso cuando se pretende realizar 
estudios de vulnerabilidad y riesgo de edificios, puentes o lineas vitales en dichos sitios. 
En tales situaciones, resulta deseable disponer de algun metodo que permita realizar la 
extrapolation de la information instrumental de manera liable. Igualmente, la estimacion 
de los contornos de igual valor de una variable macroslsmica dada resulta de importancia 
para la modelacion de la tectonica regional. 

Es bien sabido que la atenuacion de las medidas de la intensidad sismica mencionadas an- 
teriormente se describe usualmente de forma estadlstica por medio de ecuaciones de regresion 
no lineal. Por ejemplo, una funcion usada tipicamente para caracterizar la atenuacion de la 
aceleracion pico es (Newmark y Rosenblueth, 1976) 

®max = Cj exp (c 2 M)(R + c 3 ) c 4 (1) 

donde M es la magnitud, R es la distancia epicentral y c { ,i = 1, . . . ,4 son los coeficientes 
de regresion. A pesar de la gran utilidad mostrada por ecuaciones de este tipo para el 
analisis de la amenaza (tambien llamada peligrosidad) sismica, resultan de poco valor para 
estimar la aceleracion maxima ocurrida en un sitio carente de registro instrumental. Esto 
se debe a que la forma de la ecuacion es la de una elipse perfecta, la cual tiene siempre 
sus ejes orientados en las direcciones horizontal y vertical. Esto contrasta con la forma 
irregular de los contornos observados de iso-aceleracion, los cuales, ademas, presentan unos 
ejes fuertemente determinados por la orientacion de la zona de ruptura. 

Para la estimacion de la atenuacion espacial de cada evento se requiere, en consecuencia, 
de metodos diferentes. En principio, cabe hacer uso de las tecnicas establecidas de estadlstica 
espacial (Ripley, 1981; Cressie, 1995), las cuales han sido desarrolladas por algunas ciencias 
de la Tierra, tales como la mineria y la agricultura y han sido luego aplicadas en disciplinas 
de lo mas diverso, tales como la demografia, la ciencia de los materiales, la sociologia, etc. En 
esencia, estas tecnicas apuntan a la obtencion de la estructura de la covarianza de los datos 
espaciales, sobre cuya base resulta posible estimar los valores de la variable bajo estudio 
en los sitios carentes de medicion. En vista de que la covarianza esta ligada por medio de 
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una transformacion bidimensional de Fourier a la densidad espectral de potencia, resulta 
posible ademas utilizar tecnicas de simulacion para estimar los datos ausentes de manera 
probabilista. 

En este articulo se explora una alternativa diferente para este proposito, consistente en la 
aplicacion de un tipo de dispositivo de inteligencia artifical, conocido como redes neuronales 
multi-estrato, para la identification de las tendencias de atenuacion espacial de los datos 
disponibles, a partir de lo cual se podria inferir los valores de la variable en otros sitios. 
Este uso de las redes neuronales se encuentra justificado en el hecho de que estas se han 
convertido en una tecnica estandar en el area del Reconocimiento de Patrones (Ripley, 1996; 
Theodoridis y Koutroumbas, 1999), la cual tiene por objetivo central la clasificacion de datos 
dispersos, que es justamente lo que se pretende con el calculo de contornos espaciales. En 
general, las redes neuronales hacen parte de la Inteligencia Artificial (Jang et ai, 1997)- un 
area de rapido desarrollo que comienza a encontrar aplicacion en multiples campos, entre 
los cuales se encuentran la Sismologia y la Ingenieria Sismica (Garcia y Zonno, 1996). Las 
principals ventajas de este enfoque son las siguientes: 

1. No se requiere hacer hipotesis sobre la estructura probabilista de los datos. En 
contraste, otras tecnicas de Reconocimiento de Patrones tales como la clasificacion Bayesiana 
requieren suposiciones de este tipo (Ripley, 1996), al igual que en la Estadistica Espacial se 
requiere suponer que los datos responden a una distribution de Gauss (Cressie, 1995). 

2. Se pueden analizar patrones complejos sin dificultad evidente, ya que el perceptor neu- 
ronal presenta propiedades de asociacion paralela de rasgos que no requieren compromises 
entre los datos y alguna funcion prescrita, como sucede en las tecnicas de regresion. _ 

La ultima caracteristica de las mencionadas explica el exito obtenido por Emami et al. 
(1996) en el calculo de la atenuacion de la aceleracion maxima como funcion de la magnitud 
y la distancia epicentral, como alternativa al uso de ecuaciones del tipo (1). Tal aplicacion 
es una buena muestra del uso potencial de estas tecnicas en estudios globales de amenaza 
sismica. En el caso presente se busca examinar la posibilidad de emplear las redes neuronales 
para inferir el patron de atenuacion espacial de un suceso sismico dado sobre la base de la 
information disponible, ya que la capacidad de las redes neuronales para tratar con sistemas 
altamente no lineales las hace ideales para una situation tan compleja como la que se da en la 
dispersion de las ondas sismicas a traves de la corteza terrestre. Para este estudio se entreno 
una red neuronal con los datos de aceleracion maxima registrada de las tres componentes 
de un sismo ocurrido en Colombia en 1995, buscando un error medio menor del 6 por ciento 
entre las aceleraciones registradas y las estimadas por la red. 

El articulo esta dividido en dos secciones. En la primera, se hace una presentation sucint, a 
de la arquitectura de la red neuronal multi-estrato asi como una sintesis del algoritmo de 
entrenamiento. En la segunda se presenta el analisis del caso mencionado, asi como las 
comparaciones realizadas con una medida indirecta de la dispersion espacial de la energia, 
cual es el mapa de isositas (contornos de igual intensidad de Mercalli). 
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2. CONCEPTOS BASICOS SOBRE REDES NEURONALES 

Una neurona biologica real esta compuesta de una celda y varias dendritas cjue la unen 
con otras neuronas. La actividad de cada una abarca las siguientes fases: (a) la reception 
de impulsos electricos enviados a ella; (b) La transformation de los impnlsos recibidos y (c) 
el envlo de esta informacion transformada a otras neuronas. En la imitation algoritmica 
de este proceso empleada en el campo de la Inteligencia Artificial, los pasos anteriores se 
convierten en los siguientes: 

a) La informacion enviada por la neurona 1 a la neurona m, x, , se multiplica por 
un factor de peso w ljn que representa la importaneia de aquella en la actividad de esta. A 
continuation se suma la informacion total que llega a la neurona m desde todas las neuronas 
conectadas con ella: 



L 

ym=Y, W lm X lm ( 2 ) 

1=1 

b) Se transforma el resultado por medio de una funcion no lineal 

y 3 = f(yj) 

Para este proposito la funcion mas usada es, quizas, la llamada funcion sigmoide unipolar , 
dada por 



/(*) 



l_ 

a + e -71 



( 3 ) 



c) La informacion transformada se envia a otras neuronas k, afectando cada recorrido 
por el factor de peso correspondiente w jk . 




Figura 1 Representacion esquematica de una red neuronal 

De esta manera resulta posible const.ruir un perceptor artificial completo al agrupar las 
neuronas en estratos, de la manera que muestra la figura 1, la cual corresponde a un tipo 
de red de arnplia aplicacion conocido como perceptor de tres estratos. Este dispositivo se 
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utiliza en general para aprender las relaciones de entrada-salida de sistemas complejos de 
diversas clases, tales como las examinadas por Yagawa y Okuda (1997) en el campo de la 
Mecanica de Estructuras. El proceso de entrenamiento es el siguiente: Se presenta a la red 
la informacion de un conjunto de datos a traves de la primera capa en forma normalizada. 
En este paso, los valores de las conexiones w tj , que son el objetivo del entrenamiento, 
toman usualmente valores aleatorios. Despues de sumar la informacion recibida y aplicar 
la transformacion no lineal en cada neurona de la capa intermedia (u oculta) se determinan 
los valores de salida en la capa del mismo nombre, aplicando las mismas operaciones de 
suma y transformacion. Estos resultados se comparan con los reales, con los cuales hay, 
necesariamente, una diferencia, la cual se utiliza para actualizar los valores de las conexiones. 
El proceso se repite para el siguiente conjunto de datos de entrada-salida hasta agotar la 
base de datos disponible. Lo normal es que al final de este proceso, llamado epoca, el error 
no se encuentre dentro de lo deseable, por lo cual se deba repetir de nuevo el proceso para 
multiples epocas hasta que se alcance dicho objetivo. 

Este algoritmo se conoce con el nombre de propagation hacia atras. Los siguientes son 
los pasos mas relevantes del entrenamiento de acuerdo con uno de los metodos mas usados, 
llamado la regia delta (Cichocki y Unbehauen, 1993): 

1. Se aplica un vector de datos \x x ,x 2 , . . . ,x L ] en las L neuronas del primer estrato. 
Usu alm ente se introduce una neurona adicional en esta capa, llamada de sesgo, cuyo valor 
tipico es x L+l = 1. 

2. Se calcula la informacion que llega a todas las neuronas del segundo estrato de acuerdo 
con la ecuacion (2). En particular, 



L 

Vm = Y. w im x t +^ m 

(=i 

donde w m es el peso de la neurona de sesgo. 

3. Se aplica la transformacion nolineal a la entrada total de cada neurona m del segundo 
estrato, m — 1,2 M: 



llm=fm{ym ) ( 5 ) 

En esta ecuacion el sub-indice de la funcion insinua la posibilidad de tener diferentes 
funciones de transformacion para cada neurona, lo cual, empero, no es el caso corriente. 

4. Se calcula la entrada y la salida de la tercera capa de manera semejante a lo hecho en 
la segunda: 



M 



Z n=Y. W mnV m + W n 

m— 1 



( 6 ) 



= /nOO 



( 7 ) 



5. El error de estimation que media entre el dato real de cada neurona n, r n , y el 
estimado por la red es 



— r — 



( 8 ) 



En la regia delta se busca minimizar el error total cuadratico de la capa de salida, dado por 



54 



J. E. Hurtado y M. Builes 



e 



N 



T8 2 

/ -j n 



n— 1 



( 9 ) 



con respecto a los pesos w mn . Esto implica calcular la derivada parcial de la funcion de error 
con respecto a w mn , lo que da 



de 



= “( r n ~K) 



dfM dz n 



( 10 ) 



V " dZ n d W _ 

En el caso de la sigmoide unipolar, la derivada parcial contenida en el lado derecho de la 
ecuacion anterior es 



9f(x) = d f(x) 
dx Ax 



= 7/(*)(l - /(*)) 



Esto implica que 



JK 

dw m 



dfn(Zn) 

dz n 
d 



= 7*»(1 -*„)> 



M 



dw„ 



W mny m +™ n = Vn 



mn m= 1 

Por tanto, el gradiente de la funcion de error es 



de 

i 



= -(r„ -2 n )7i n (l ~z n )y 



(ID 

( 12 ) 

( 13 ) 

( 13 ) 



Como el error minimo se encuentra en la direccion negativa del gradiente. los pesos de 
la capa de salida en el paso [k + 1] de entrenamiento se calculan anadiendo una porcion del 
gradiente negativo a los pesos calculados en el paso k: 



w mn[k + 1] = w mn [k] + T](r n -z n )yz n (l - z n )y m (14) 

En esta ecuacion y es el llamado factor de velocidad de aprendizaje. cuyo valor se escoge 
usualmente en el rango (0, 1). Se puede agregar un factor adicional a. llamado de momentum 
para acelerar el aprendizaje (Cichocki y Unbehauen, 1993) . 

6. La actualizacion de los pesos de la capa oculta se efectua de manera similar: 

N 

W lmi k + !] "“taW+TOj 1 -Vm) X l ~ KhKiX ~ ^J W mn' k ] ( L 5) 

71=1 

7. Se repite el proceso para nuevos valores de entrada-salida (x r r n ). I = 1.2. ... L. n = 
1,2 ,N durante varias epocas hasta que el error e este dentro de lo admisible. 

8. Una vez entrenada la red se puede usar para estimar la salida correspondiente a 
entradas no presentadas previamente por medio de las ecuaciones (4) a (7). 

Existen muchas variantes del algoritmo anterior asi como varias propuestas de funciones 
de activacion, las cuales pueden ser ensayadas para mejorar la convergencia. En el presente 
caso, sin embargo, el algoritmo estandar descrito mostro desde el comienzo una tendencia 
uniforme hacia la convergencia en los tres casos analizados, por lo que no se considero 
necesario probar otras variantes. 
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3. MODELACION NEURONAL DE LA ATENUACION SISMICA 

Para este trabajo se aplico el algoritmo anterior a los datos instrumentales disponibles 
de un evento sismico de Colombia, cual es el sismo de Tauramena, que tuvo lugar el 19 de 
enero de 1995. Su magnitud es 6.5, la localization del epicentro es 5.01 N, 72.95 W y la 
profundidad focal estimada es 25 km. Como ha quedado dicho, el proposito es indagar la 
aplicabilidad de la red neuronal descrita para simular los contornos de iso-aceleracion en la 
region de infiuencia del suceso sismico. La red neuronal se entreno con la informacion sobre 
longitud y latitud geograficas de cada estacion de registro en las nueronas de entrada y las 
aceleraciones maximas, en valor absolute, de las tres componentes (EW, vertical y NS), en 
las de salida. Se hizo uso de la neurona de sesgo en la capa de entrada y se normalizaron 
los datos, tanto de entrada como de salida en el rango [0, 1]. La comparacion con los datos 
reales registrados que se muestra mas adelante corresponde, obviamente, a las estimaciones 
de la red una vez retirada la normalization. 

La informacion instrumental disponible en este caso (INGEOMINAS, 1995a, 1995b) 
comprende 16 registros. Las aceleraciones maximas horizontales, independientemente de 
la direction, aparecen en la figura 2. Un examen superficial de esta figura permite concluir 
que resulta dificil trazar los contornos buscados debido a cierta naturaleza contradictoria 
de los datos. La figura 3 muestra las historias de la media aritmetica y de la media mas 
una desviacion estandar del error de estimacion de la red. En total se necesito un mimero 
aproximado de 1,500,000 epocas para lograr un error medio inferior al 5.5 %. La figura 4 
muestra los contornos de iso-aceleracion obtenidos usando la red ya entrenada con datos 
de entrada de latitud y longitud correspondientes a una malla de puntos geograficos de 
ancho 0.1 grados. La comparacion con la figura 2, en la que los datos instrumentales 
se encuentran en una franja de tendencia horizontal, parece insinuar que los contornos 
predichos, de tendencia claramente vertical, resultan mas de la imagination de la red que 
de los datos disponibles. Para poder aseverar la bondad de la estimacion se requiere, en 
consecuencia, usar los contornos de Intensidad de Mercalli, que se calculan corrientemente 
en estos casos, debido a la relation directa (aunque no proporcional) que se presume entre la 
intensidad de Mercalli y la aceleracion maxima horizontal. En consecuencia, los contornos 
de aceleracion deberian mostrar en principio cierto paralelismo con los de intensidad. Para 
el sismo en mention, dichos contonos (Escallon y Martinez, 1997) aparecen en la figura 5 
sobre una ventana geografica igual a la usada en la figura 4. (En este respecto cabe observar 
que, debido a la incapacidad reconocida de las redes neuronales de extrapolar fuera del 
rango de los datos de entrenamiento (Cichocki y Unbehauen, 1993), los contornos han sido 
calculados sobre una ventana menor que la de la figura 2). La orientacion de los contornos 
de isosistas mostrados en la figura 5 corrobora la orientacion vertical de los contornos de 
atenuacion estimada por la red. La semejanza de los contornos de igual intensidad con los 
de iso-aceleracion estimados por la red resulta evidente. 

Para completar la informacion de los calculos realizados con la red, la figura 6 muestra 
los contornos de igual aceleracion vertical del sismo. Puede verse que todos siguen de igual 
forma la tendencia vertical mencionada. 




Longitud 
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Figura 4 Estimation de la aceleracion maxima horizontal dada 
por la red neuronal 



Intensidades de Mercalli - Sismo de Tauramena 




Figura 5 Isosistas del sismo estudiado (Escallon y Martinez 1997) 
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Figura 6 Estimacion de la aceleracion maxima vertical dada por 
la red neuronal 



CONCLUSIONES 

En este artlculo se ha explorado la posibilidad de usar las redes neuronales conocidas como 
perceptores multi estrato para estimar las curvas de igual aceleracion maxima producida 
por movimientos sismicos. Para ello se han usado los datos correspondientes a un sismo 
ocurrido en Colombia del cual se tiene la informacion de contornos regionales de igual 
intensidad de Mercalli ademas de buena informacion instrumental. Se ha mostrado que todos 
los contornos de maxima aceleracion estimados por la red (dos componentes horizontales 
y la vertical) tienen tendencias semejantes a los de las isosistas. Por el contrario, la 
estimacion de los mismos a partir de ecuaciones estadisticas de atenuacion esta rodeada 
de incertidumbre en este respecto, debido a la rigidez de las ecuaciones (es decir, a su forma 
matematica preestablecida) y a la independencia que imponen a la estimacion con respecto 
a las coordenadas geograficas del sitio, lo cual es contrario a la evidencia. 
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RESUMEN 

Se analizan las zonificaciones sismicas y los espectros de diseno elastico, de las normativas 
vigentes al ano 2000, en Venezuela, Colombia, Ecuador y Peru. Previamente se comparan 
las zonificaciones obtenidas en los indicados paises, a principios de la decada de los ochenta, 
cuando se empezo a calcular la peligrosidad sismica en terminos probabilisticos, con las 
actuates 

Los mapas de zonificacion que se analizan determinan la aceleracion maxima del suelo en 
roca, para una vida util de la estructura de 50 anos y con una probabilidad de excedencia 
del 10% y son la base para la obtencion de los espectros elasticos de diseno. 

Se destaca la necesidad de contar con un nuevo formato de presentacion de la zonificacion 
sismica, de tal manera que permita encontrar la aceleracion maxima del suelo para diferentes 
periodos de retorno de la estructura. 



ABSTRACT 

The seismic zonification and elastic design spectra of the current regulations, updated 
to the year 2000, in Venezuela, Colombia, Ecuador and Peru are analyzed. Previously, the 
zonifications obtained in the same countries, using the regulations of the early eighties, when 
seismic danger started being calculated in probabilistic terms, were compared to current 
regulations. 

The zonification maps analyzed present the maximum soil acceleration in the rock, tor a 
useful structure life of 50 years, and an excedence probability of 10%. This is the base to 
obtain elastic design spectra. 

The need to have a new presentation format of the seismic zonification is remarked, m 
such a way to allow finding the maximum soil acceleration for different return periods of 
the structure. 
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1 INTRODUCTION 

Existen algunas formas de cuantificar y presentar los mapas de zonificacion sismica de 
una determinada region. Una de ellas es en base al analisis de las curvas isosistas de 
eventos regist.rados en el pasado (15) . Para ello, los danos o efectos, producidos por un sismo 
son evaluados en funcion de la Intensidad. Con esta variable, como base se han realizado 
algunas zonificaciones a nivel de Pais, como es el caso de la primera norma sismica de Peru 
de 1971 o a nivel regional^ 12 *. 

Otra forma, de zonificar un Pais, aunque un tanto mas elaborada es por medio de los 
mapas sismotectonicos, que incluyen: information tectonica, geologica, geofisica, geotecnica 
v sismica 121 b Desde el punto de vista de Ingenieria Sismo resistente, esta zonificacion no 
aporta con parametros con los cuales se pueda determinar las fuerzas de diseno sismico de la 
estructura. Por esta razon, la corriente mas utilizada para zonificar sismicamente un Pais, 
es la que se indica a continuation. 

En la tercera forma de calculo, se procede a evaluar la peligrosidad sismica en terminos 
probabilisticos (6 ' 8,16) . Para la evaluation, es necesario definir leyes de atenuacion para los 
movimientos del suelo. El estudio sismotectonico indicado en el parrafo anterior sirve para 
determinar las areas fuentes o zonas sismogeneticas. En esta corriente de evaluacion, el ATC 
3-06 plantea la construction de espectros de aceleracion elastica en base a la aceleracion 
maxima del suelo en roca, la misma que se obtiene de mapas de regionalization sismica 
originados en curvas de igual aceleracion para una vida util de la estructura de 50 anos y 
con una probabilidad de excedencia del 10%. 

Las zonificaciones sismicas que se han realizado en Venezuela, Colombia, Ecuador y Peru, 
teniendo en cuenta la vida util y probabilidad de excedencia indicados se analizan en los 
siguientes apartados. 



2 ZONIFICACION SISMICA DE VENEZUELA 

En la figura 1, se indica lo que podria ser, el primer mapa de zonificacion sismica 
de Venezuela, que fue publicado en la Norma COVENIN 1756-82 (13) . Este trabajo fue 
realizado en el Instituto de Materiales y Modelos Estructurales por Jose Grases* 20 \ para el 
efecto se implemento en un computador H.P. 1000 E., el programa “Seismic Risk Analysis”, 
desarrollado por Robin McGuire (22) en 1975. 

Se destaca que se han definido cuatro zonas sismicas, cuya aceleracion maxima del suelo 
en roca, varia desde 0.08 g (8% de la aceleracion de la gravedad) para la zona uno, hasta 
0.30 g, para la zona cuatro. Las zonas dos y tres estan caracterizadas por aceleraciones de 
0.15 g. y 0.22 g, respectivamente. 

En 1998, se publica la Norma Covenin 1756-98 (14) y en ella se presenta la nueva zonifi- 
cacion sismica de Venezuela, que podria considerarse como la ultima publicada a nivel de 
normativa para edificaciones sismorresistentes, la misma que se indica en la figura 2. 

En la Norma Venezolana de 1998, se clasifica el Pais, con fines de Ingenieria Sismica, 
en siete zonas, que van desde 0.10 g, en la zona uno hasta 0.40 g, en la zona de mayor 
peligrosidad sismica, que es la siete. 

Los cambios mas significativos que se han dado en los dos mapas, indicados en las figuras 
1 y 2, son que en la Norma Covenin 1756-98 la zona 1, esta caracterizada por una aceleracion 
de 0.10 g, antes en 1982, la aceleracion era nula. Por otra parte, en las zonas seis y siete de 
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Figura 1 Zonificacion sfsmicade Venezuela, segun COVENIN 1756-82. 




Figura 2 Zonificacion sfsmicade Venezuela, segiin COVENIN 1756-98. 



la Norma de 1998, se espera una aceleracion maxima del suelo en roca de 0.35 g. y 0.40 g, 
en la norma de 1982 estas dos zonas estaban definidas por una peligrosidad sfsmica de 0.3g. 

Otra diferencia que existe entre las dos normas de 1982 y 1998, en cuanto se refiere 
a la determinacion de la aceleracion maxima del suelo en roca, es la incorporacion en la 
Norma Covenin 1756-98 del factor de correction <p, el cual modifica la aceleracion, indicada 
en la figura 2, en funcion del tipo de suelo. En efecto, en suelos en los cuales hay menos 
probabilidad de amplification de las ondas sfsmicas, como se tiene en los suelos duros, se 
puede disminuir hasta en un 15%, la aceleracion especificada en la figura 2. 
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Figura 3 Zonificacion si'smica de Colombia, segun ACIS 100-81. 



3 ZONIFICACION SfSMICA DE COLOMBIA 

Cuando se habla de Peligrosidad Sismica de Venezuela, es importante estudiar los tra- 
bajos desarrollados por Jose Grases^ 1,19 *, de igual forma cuando este tema es tratado sobre 
Colombia, se deben conocer las contribuciones realizadas por Alberto Sarria (26 h 

En efecto, en Colombia Alberto Sarria ha sido autor o coautor de destacados trabajos, 
uno de ellos fue el que dirigio a Jose Atuesta' 10 ' en la Universidad de los Andes, como 
proyecto de grado, en 1972. Del resultado de esta investigacion Colombia se clasifica en 
cuatro zonas definidas de cero a tres. A1 analizar este trabajo con la Normativa Sismo 
Resistente NSR-98 (24 1, se observa una buena similitud, la zona cero corresponde a la que en 
la normativa NSR-98 la denominan “baja” y la zona 3 a la que denominan “alta”. 

Antes de la aparicion del primer codigo de construcciones sismorresistentes, en 1984, la 
Asociacion Colombiana de Ingenieria Sismica, publico en 1981 la Norma ACIS 100-81 (2) la 
misma que fue utilizada de manera voluntaria por una gran cantidad de profesionales de la 
ingenieria. En la figura 3, se indica la zonificacion presentada en la Norma ACIS 100-81, 
se observa que al Pais se lo ha dividido en 6 zonas, en la zona 1, la aceleracion maxima en 
roca es nula y en la zona 6, la aceleracion tiene un valor de 0.30 g. 





